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Dimers and cluster integrable systems 

A.B. Goncharov, R. Kenyon 



Abstract 



We show that the dimer model on a bipartite graph T on a torus gives rise to a quantum 
integrable system of special type, which we call a cluster integrable system. The phase 
space of the classical system contains, as an open dense subset, the moduli space £r of 
line bundles with connections on the graph I\ The sum of Hamiltonians is essentially the 
partition function of the dimer model. 

We say that two such graphs Ti and T 2 are equivalent if the Newton polygons of the 
corresponding partition functions coincide up to translation. We define elementary transfor- 
mations of bipartite surface graphs, and show that every two equivalent minimal bipartite 
graphs are related by a sequence of elementary transformations. For each elementary trans- 
formation we define a birational Poisson isomorphism C? 1 — > Cr 2 providing an equivalence 
of the integrable systems. We show that it is a cluster Poisson transformation, as defined in 

We show that for any convex integral polygon N there is a non-empty finite set of minimal 
graphs r for which N is the Newton polygon of the partition function related to I\ Gluing 
^ ■ the varieties Cr for graphs T related by elementary transformations via the corresponding 

cluster Poisson transformations, we get a Poisson space Xpf. It is a natural phase space for 
the integrable system. The Hamiltonians are functions on Xpj, parametrized by the interior 
points of the Newton polygon N . We construct Casimir functions whose level sets are the 
symplectic leaves of Xjy ■ 

The space Xm has a structure of a cluster Poisson variety. Therefore the algebra of regular 
functions on Xm has a non-commutative q-deformation to a *-algebra O q (Xn). We show 
that the Hamiltonians give rise to a commuting family of quantum Hamiltonians. Together 
with the quantum Casimirs, they provide a quantum integrable system. Applying the general 
quantization scheme [FG2j . we get a ^-representation of the *-algebra O q {X^[) in a Hilbert 
space. The quantum Hamiltonians act by commuting unbounded selfadjoint operators. 

For square grid bipartite graphs on a torus we get discrete quantum integrable systems, 
where the evolution is a cluster automorphism of the *-algebra Oq(Xjj) commuting with the 
quantum Hamiltonians. 

Any graph G on a torus T gives rise to a bipartite graph Tq on T. We show that the phase 
space X related to the graph Tq has a Lagrangian subvariety H, defined in each coordinate 
system by a system of monomial equations. We identify it with the space parametrizing 
resistor networks on G. The pair (X, TV) has a large group of cluster automorphisms. In 
particular, for a hexagonal grid graph we get a discrete quantum integrable system on X 
whose restriction to TZ is essentially given by the cube recurrence of [GSj . 

The set of positive real points X^ (K>o) of the phase space is well defined. It is isomorphic 
to the moduli space of Harnack curves with divisors studied in |KO] . The Liouville tori of 
the real integrable system are given by the product of ovals of the Harnack curves. 

In the sequel |GK] to this paper we show that the set of complex points Xn (C) of the 
phase space is birationally isomorphic to a finite cover of the Beauville complex algebraic 
integrable system related to the toric surface assigned to the polygon N. 
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1 Introduction 

A line bundle V on a graph T is given by assigning a 1-dimensional complex vector space V v to 
each vertex of T. A connection on V is a choice of isomorphism <p vv / : V v — > V v > whenever v, v' 
are adjacent, satisfying <j> v > v = <f>~ v i- 

We construct a Poisson structure on the space of line bundles with connections on a bipartite 
graph embedded on a surface. When the graph is embedded on a torus T we construct an 
algebraic integrable system, with commuting Hamiltonians which can be written as Laurent 
polynomials in natural coordinates on Cy , the moduli space of line bundles with connections on 
r. The Hamiltonians are sums of dimer covers of V. 

A dimer cover of a graph T is a set of edges with the property that every vertex is the endpoint 
of a unique edge in the cover. Probability measures on dimer covers have been the subject 
of much recent work in statistical mechanics. In this paper we study some non-probabilistic 
applications of the dimer model, in particular to the construction of cluster integrable systems. 

Below we introduce the key non-technical definitions which we use throughout the paper, 
and discuss the main results. 

1.1 Dimer models and Poisson geometry. 

A surface graph V is a graph embedded on a compact oriented surface S whose faces, i.e. the 
connected components of S — V, are contractible. It is homotopy equivalent to an open surface 
So C S, obtained by putting a puncture in every face. A bipartite graph is a graph with vertices 
of two types, black and white, such that each edge has one black and one white vertex. 

1.1.1 Conjugate surface graph and Poisson structure 

Let Cr be the moduli space of line bundles with connections on a graph T. Any oriented loop L 
on r gives rise to a function Wl on £r given by the monodromy of a line bundle with connection 
along L. One has W-l = W7 , where —L is the loop L with the opposite orientation. The 
monodromies around the loops provide an isomorphism 

/: r ^Hom(Fi(r,z),c*) = J H' 1 (r,c*). (i) 

For any surface graph T, the moduli space Cr has the standard Poisson structure. Namely, 
since T is homotopy equivalent to the surface Sq, we have Cr = ^(Sq^C*). This latter space 
has a Poisson structure provided by the intersection pairing on Hi(So,7*). 

Now let r be a bipartite surface graph. Then we introduce a new Poisson structure on £p 
related to the standard one by a "twist". 

A surface graph is the same thing as a ribbon graph: a ribbon graph is a graph with an 
additional structure given by, for each vertex, a cyclic order of the edges at that vertex. Given 
a ribbon graph, we can replace its edges by ribbons, getting an oriented surface. 

A bipartite surface graph T gives rise to a new ribbon graph r^, = T, the conjugate graph, 
obtained by reversing the cyclic orders at all white vertices. The topological surface with 
boundary corresponding to the ribbon graph T w is called the conjugated surface S w . 



1 Reversing the cyclic order at the black vertices we get a ribbon graph Fb, which differs from F w by the 
orientation. 



Evidently a line bundle with connection on the graph T is the same thing as a line bundle 
with connection on the original graph T. So there is a canonical isomorphism 

C T = C f . (2) 

The standard Poisson structure on the moduli space £p combined with the isomorphism ([2]) 
provides the Poisson structure on Cr we need. Precisely, consider the intersection pairing on 

e §w :H 1 (S w ,Z)AH 1 (S w ,Z)^Z. (3) 

Given two loops Li,L 2 on T, we define the Poisson bracket {Wl 1 ,Wl 2 } by setting 

{W Ll ,W L2 } = e Sw {L 1 ,L 2 )W Ll W L2 . (4) 

and extending it via the Leibniz rule to a Poisson bracket on the algebra generated by the 
functions Wl , which coincides with the algebra of regular functions on the variety £r ■ 

Natural coordinates on C-p- The orientation of the surface S induces orientations of the 
faces of r. The group Hi(T, Z) is generated by the oriented boundaries d(F) of the faces F of T, 
whose orientation is induced by the orientation of F, and loops generating Hi(S,Z). The only 
relation is that the sum of boundaries of all faces is zero. There is an exact sequence, where the 
map 7r is induced by the embedding r *->■ S: 

— >Kerir — ►fl , 1 (T;Z) -^ #i(5;Z) — >■ 0. 

The subgroup Ker tt is generated by the oriented boundaries d(F) of faces F. 

Summarising, the monodromies of connections around all but one oriented faces F of T, 
augmented by the monodromies around loops generating ^(S 1 , Z), provide a coordinate system 
{Xi} on the moduli space Cr- The Poisson structure in these coordinates has a standard 
quadratic form 

{Xi,Xj} = EijXiXj, Eij € Z, (5) 

where e is the intersection pairing on S w . 

Zig-zag loops and the center of Poisson algebra of functions on Cy- A zig-zag path on 
a ribbon graph graph T is an oriented path on V which turns maximally left at black vertices 
and turns maximally right at black vertices; see Figure [TJ There is a bijection 

{zig-zag loops on T} <-> {boundaries of the holes on S w }. (6) 

The orientations of zig-zag loops match the orientations of the boundary loops on S w induced 
by the orientation of S w . 




Figure 1: A Zig-zag path on a ribbon graph. 



Lemma 1.1. Let Z be an (oriented) zig-zag loop on a bipartite oriented surface graph T. Then 
as Z runs over zig-zag loops, the functions Wz generate the center of the Poisson algebra 0{Cy) 
of functions on Cr- The product of all of them is 1. This is the only relation between them. 

We call the functions Wz the Casimirs. 

Proof. It follows from the canonical bijection ([6]). Evidently the boundary loops on S w generate 
the kernel of the intersection pairing on H\{S w ,1j). Clearly the sum of them is zero, and this is 
the only relation between them. □ 

There is a similar Poisson moduli space £r for a bipartite graph T on an oriented surface S 
with boundary. 

1.1.2 Gluing the Poisson phase space 

The set of homology classes of dimer coverings of T (as defined below) defines a convex polyhe- 
dron iV C Hi(S, M). It is the Newton polyhedron of the dimer partition function. Two bipartite 
surface graphs are equivalent if their Newton polyhedra coincide up to translation. In each 
equivalence class there are minimal bipartite surface graphs, which have the minimal number of 
faces. We prove that every two minimal bipartite surface graphs are related by certain elemen- 
tary transformations, consisting of "spider" moves and "shrinking of a 2-valent vertex" moves. 
Each elementary transformation T\ — > T2 gives rise to a cluster Poisson birational isomorphism 

ir u r 2 : C Tl -^ £r 2 (7) 

We define a Poisson variety Xn, the phase space, by gluing the moduli spaces £r via transfor- 
mations (J7j). The variety Xn depends only on N. 

Since the sets of zig-zag loops on equivalent bipartite surface graphs are naturally identified, 
the gluing maps (J7|) preserve the Casimirs. Therefore the center of the Poisson algebra of regular 
functions on the phase space Xn is also described by the Casimirs. 

1.2 Cluster integrable systems for dimer models on a torus. 

1.2.1 Classical integrable system 

Let us restrict to the case when T is a bipartite graph on a torus T. In this case the Newton 
polyhedron N is a convex polygon in the plane Hi(T,R) with vertices at the lattice H\(T, Z). 
Zig-zag loops are parametrized by the intervals on the boundary dN between two consecutive 
integral points, called primitive boundary intervals. 

We show that the partition function of the dimer model on T gives rise to a regular function 
V on Xn, the modified partition function, well defined up to a multiplication by a monomial of 
Casimirs. It has a natural decomposition into a sum of components H a , called the Hamiltonians, 
parametrized by the internal integral points of the Newton polygon N. These Hamiltonians are 
Laurent polynomials in the X^s. The Hamiltonian flows are well defined. 

Theorem 1.2. The Hamiltonian flows of the H a commute, providing an integrable system on 
Xn . Precisely, we get integrable systems on the generic symplectic leaves of Xn, given by the 
level sets of the Casimirs. 



In particular, 

dimAjv = 2t(JV) + e(N) - 1, 

where i(N) (respectively e(N)) is the number of integral points inside (respectively on the 
boundary) of the polygon N. 

1.2.2 Quantum integrable system 

The algebra of regular functions on any cluster Poisson variety X has a natural (/-deformation 
to a *-algebra O q (X), defined in Section 3 of |FG1| . It can be realized, in many different ways, 
as a subalgebra in the quantum torus algebra. The latter is generated by the elements X^ 
satisfying the following relations, where q is a formal variable commuting with Xj's: 

q-^XiXj = q-^XjXi, e tj € Z. (8) 

It has an involutive antiautomorphism *, acting on the generators by 

*{Xi)=Xi, *{q) = q-\ 

So in our case we get a *-algebra O q {Xn)- The s^ in (jHJ) is the same as the Poisson tensor 
in (J3J) . The quantum analog of Theorem 11.21 is 

Theorem 1.3. The Hamiltonians H a give rise to commuting selfadjoint elements M a in the *- 
algebra O q (X^), which we call quantum Hamiltonians. The Casimirs give rise to the generators 
of the center of the algebra. 

Let us assume now that \q\ = 1. Then, according to [FG2], for any cluster Poisson variety 
X, there is a family of Hilbert spaces TL X , parametrized by the characters x °f the center of the 
algebra O q (X), and a ^representation of the *-algebra O q (X) in (a certain Schwartz space S x , 
dense in) M x . 

So in our case we get a family of Hilbert spaces TLn, x -, parametrized by a collection of 
real numbers x> the eigenvalues of the Casimirs. The quantum Hamiltonians are realized by 
commuting unbounded selfadjoint operators in Hn,x- 

Problem 1. Find the generalized eigenf unctions of the quantum Hamiltonians, and the spectral 
decomposition. 

We call the emerging object a cluster integrable system. A key point is that a cluster inte- 
grable system admits a quantum version, as described above. We believe that many interesting 
integrable systems are cluster integrable systems. 

1.2.3 Discrete cluster integrable systems. 

In the case of the dimer model on a square grid on a toruaj our quantum integrable system 
has a remarkable automorphism A, called the octahedron recurrence or, in other language, 
Hirota's bilinear difference equation (HBDE) [HirJ. We show that it commutes with the quantum 
Casimirs and quantum Hamiltonians, and is given by an element of the cluster modular group 
[FG1] of the cluster variety X - see Section El 



2 The square grid dimer models on a torus are parametrized by sublattices of the isomorphism group 1? of the 
standard square grid bipartite graph in R 2 . 



There are similar discrete cluster integrable systems related to hexagonal graphs G on a 
torus, where the discrete evolution operator is described by the cube recurrence, see Sections 
I6.2H6.3I Its phase space is the phase space of the dimer model related to the resistor network 
model on G, see Section 11.31 

We stress that discrete cluster integrable systems are discrete quantum integrable systems. 

1.3 Cluster nature of a resistor network model on a torus 

1.3.1 Dimer model arising from a graph on a torus 

A special case of the dimer model is one arising from a planar resistor network G via the 
generalized Temperley's trick [KPW]. Namely, an arbitrary graph G on a torus gives rise to 
a bipartite graph T = Tq, obtained by taking a barycentric subdivision of the graph G and 
declaring the centers of the edges of G to be the white vertices and the centers of the faces of G 
together with the vertices of G to be the black vertices. The graph T contains the dual graph 
G' to the graph of G, see Figure El where G is shown by a solid lines, and G' by dotted lines. 



Figure 2: The bipartite graph T^ (dashed and solid edges) assigned to a graph G (solid edges). 

Given a surface graph G embedded on a surface S there is a graph G embedded on the 
universal cover S, which is the natural lift of G. A surface graph G embedded on surface S is 
called minimal if zig-zag paths (defined as unoriented paths which turn alternately maximally 
left and maximally right at consecutive vertices) on G do not have selfcrossings and different 
zig-zag paths do not intersect in more than one edge. 

The Newton polygon N(G) of the bipartite graph T^j is a centrally-symmetric integral poly- 
gon in Hi(T, Z) centered at the origin. It coincides with the Newton polygon of the spectral 
curve of the Laplacian related to a generic resistor network on G, see Section [5j 

Any such polygon arises as a Newton polygon N(G) for a non-empty finite collection of 
minimal graphs G on a torus. Any two minimal graphs Gi and Gi on a torus with the same 
Newton polygon are related by a sequence of elementary transformations, called Y-A moves, see 
Figure [21 

A Y-A move G\ — > G2 induces a transformation Y\ — > T2 of the corresponding bipartite 
graphs. We decompose it into a composition of four spider moves (Lemma l5.1ip . which provides 
us with a cluster birational isomorphism 

fiy-A ■ C-Tx — > £r 2 - (9) 



Gluing the tori £p assigned to the bipartite graphs T related to the equivalent minimal graphs 
G, we get a cluster Poisson space X . |£| 



1.3.2 Resistor network Lagrangian subvariety 

The moduli space £r G of line bundles with connections on the graph Tq has a subvariety TZq 
defined by the condition that the monodromies over any loop a in G and any loop a' in G' in the 
same homology class on T are equal. These equations are given by equating certain monomials 
in the X^s to 1. 

We prove that TZq is a Lagrangian subvariety of the Poisson variety £r G , i.e. its intersections 
with the generic symplectic leaves are Lagrangian subvarieties. We show that the space IZq 
parametrizes resistor networks on G, justifying the name resistor network Lagrangian subvariety. 

Equivalently, the conjugated surface Sr G is canonically realized as the boundary of a certain 
handlebody. 

We show that the Y-A-cluster transformation ([9]) provides a birational isomorphism of the 
Lagrangian subvarieties 

TZg 1 — > 1Zg 2 ■ 






Gluing the subvarieties TZq via these birational isomorphisms, we get a Lagrangian subspace 7Z 
of the corresponding cluster Poisson space X. It gives rise to a uniquely defined up to a constant 
functional on the Schwartz space S x , see Section 15.2.71 

1.4 Toric surfaces, spectral data and the Beauville integrable system 

Let T := Tj = (C*) 2 be the algebraic torus related to the torus T. Recall that a bipartite graph 
r on T gives rise to a Newton polygon N in H\(T, R) with vertices at Hx(T, Z). Such a Newton 
polygon N, considered up to translation, provides us a toric surface M together with a line 
bundle C. The algebraic torus T acts on the pair (J\f,C). The points of M where the action is 
not free form the divisor at infinity A/"oo of M. It is a collection of P^s which is combinatorially 
isomorphic to the polygon N. 



1.4.1 The Kasteleyn operator 

A Kasteleyn line bundle with connection k on a bipartite surface graph T is a line bundle with 
connection with monodromy (— 1) ' 2+1 around faces having I edges, and ±1 around topologically 
nontrivial loops. There are 2 2g non-isomorphic choices of k when S has genus g. Given a line 
bundle with connection V on a bipartite surface graph T, we construct a new one V8k, and 
set 

V B := (V®k) v , (10) 

i>G{black vertices of T} 



V w := {V®k) v . 

-uG{white vertices of T} 

A connection on the line bundle V ® k is encoded by an operator 



(11) 



K:V B 



V 



w> 



Gluing this space we use only the graphs T which are equivalent by the Y-A-moves. We may have other 
bipartite graphs which are still equivalent as bipartite graphs. 



the Kasteleyn operator, whose matrix elements are given by parallel transports for the connection 
along the edges of T. We let 

L = detVe <g> detV w 

be the determinant line (here detV denotes the highest exterior power of a vector space V). 
The determinant det K is an element of the determinant line. 

1.4.2 The spectral data. 

Let 

T := Hom(#i(T, Z), C) = #i(T, C*) = (C*) 2 ; 

The group T is identified with the group of complex line bundles with flat connections on the 
torus T. The embedding i of the graph T into the torus T provides a free action of the group T 
on the moduli space Cp, given by * i — > * ® i*L, where L is a flat line bundle on the torus T. 
This action preserves the Poisson structure on C-p. 

The partition function can be calculated as the determinant of the Kasteleyn operator. It is 
a section of the determinant line bundle over Cp. Let us restrict it to a T-orbit on £p. Then it 
extends to a section of the line bundle C <g> L over the toric surface M which compactifies the 
orbit. The divisor of zeros of this section is called the spectral curve C. 

The spectral curve C intersects the divisor at infinity Moo at a divisor Coo- A choice of a 
point of A/"o := M — Moo provides a coordinate on every component of the divisor Moo- The 
divisors Coo satisfy a single condition: 

The product of the coordinates of all points of Coo is equal to 1. (12) 

We call such divisors admissible divisors at infinity. There is a bijection 

v : {zig-zag paths on T} ^—t divisor at infinity Coo- (13) 

It is uniquely determined by the condition that the Casimir functions on the phase space given 
by the monodromies along zig-zag paths a are proportional to the Casimir functions given by 
the coordinates of the boundary points v{a) on the corresponding components of Moo. 

Pick a black vertex b of T. A non-zero vector in the fiber Vb provides a section Sb of the sheaf 
CokerK. The latter sheaf is a line bundle for a smooth C. The restriction of its divisor of zeros 
to Mo is a degree g effective divisor S, where g is the genus of C. 

The triple (C,S,v) is called the spectral data. It was studied in [KOJ. 

1.4.3 A cover of the Beauville integrable system. 

Let S be the moduli space parametrizing the spectral data. Let B be the moduli space of pairs 
(C,v). Forgetting S, we get a projection 

ir:S^B, (C,S,v)t—*(C,v). (14) 

We claim that it is an integrable system, which covers the Beauville integrable system - the 
latter does not take into account parametrizations (|12p . 
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Namely, let S x C S be the subspace parametrizing the spectral data (C, S, u) with a given 
divisor Coo- It carries a rational symplectic form f^s defined as follows. The canonical T-invariant 
symplectic form uj on A/o gives rise to the Beauville symplectic form on Sym 9 A/o: 



^b:=E 



9 

1=1 



Here s, : M 9 — > A/" is the projection to the i-th factor. There is a projection 

p : 5 X — ► SymWo, (C, 5, i/) .— ► 5. 

It is a Galois cover, whose Galois group is a product of symmetric groups acting by permuting 
the points of C^ at each boundary component, and hence altering the bijection (fT3|) . Indeed, 
one shows that, given a generic admissible divisor at infinity C^ C TVoo an d a generic degree g 
divisor S C A/Jj, there is a unique spectral curve C passing through S and intersecting A/^ at 
the divisor Coo. So the pull back $7 := p*^s is a rational symplectic form on S x . Let B x C £> 
be the subspace given by the condition that C n A/"oo is a given divisor. There is a fibration 

^x : S x — »• B x- ( 15 ) 

Its fibers are Lagrangian. The fiber over a point (C, z/) is identified with an open part of Pic 9 (C), 
which is a principal homogeneous space over the Jacobian of C . So (|15p and therefore (|14p are 
algebraic integrable systems. 

1.4.4 Connection with the dimer integrable system. 

A proof of Theorem 11.41 below will be given in |GK| . Combined with the results of this paper it 
allows to quantize the integrable system ([14 



Theorem 1.4. For any Newton polygon N , the spectral data provides a birational isomorphism 
S over the base B: 

x A s 

\ v/ (16) 

B 

It identifies the dimer integrable system with the integrable system j!4\ ), matching the symplectic 
leaves of X and the symplectic varieties (S x ,£l). 

The map § intertwines the action of the torus TonA 1 with the action on S provided by the 
action of T on the surface A/". 

In Section 7 we prove a weaker statement: S is a finite cover over the generic part of S. It 
implies the independence of the Hamiltonians. 

1.5 Analogies between dimers, Teichmiiller theory and cluster varieties. 

Cluster Poisson varieties provide a framework for study of both (classical and higher) Teichmiiller 
theory |FG| and theory of dimers. Here is a dictionary relating key objects in these three theories. 
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Dimer Theory 


Teichmiiller Theory 


Cluster varieties 


Convex integral polygon N 


Oriented surface S 
with n > punctures 




Minimal bipartite graphs 
on a torus 


ideal triangulations of S 


seeds 


spider moves of graphs 


flips of triangulations 


seed mutations 


Face weights 


cross-ratio coordinates 


Poisson cluster 
coordinates 


Moduli space of spectral 
data on toric surface N 


Moduli space of framed 
PGL(2) local systems 


cluster Poisson variety 


Harnack curve + divisor 


complex structure on S 




Moduli space of 
Harnack curves + divisors 


Teichmiiller space of S 


positive real points of 
cluster Poisson variety 


Tropical Harnak curve 
with divisor 


Measured lamination 




Moduli space of tropical 
Harnack curve + divisors 


space of measured 
laminations on S 


real tropical points of 
cluster Poisson variety 


Dimer integrable system 


Integrable system related to 
pants decomposition of S 




Hamiltonians 


Monodromies around loops 
of a pants decomposition 





What distinguishes these two examples - the dimer theory and the Teichmiiller theory - from 
the general theory of cluster Poisson varieties is that in each of them the set of real / tropical 
real points of the relevant cluster variety has a meaningful and non-trivial interpretation as the 
moduli space of some geometric objects. 

Here by moduli space of certain objects related to the toric surface M we mean the space 
parametrizing the orbits of the torus T acting on the objects. For example, the moduli space 
of spectral data means the space «S/T. So combining results of [GK] with Theorem 14.61 would 
imply that the latter corresponds to the hypersurface X^ in the cluster Poisson variety X' N . 



Discrete cluster integrable systems in Teichmiiller theory. Consider the moduli space 
^PGh2,s of framed PGL2-local systems on S [FGJ. If S has no punctures, it is just the moduli 
space of PGL2-local systems on S, which has the classical symplectic structure. If S is a surface 
with punctures, the moduli space XpQi, 2t s nas a cluster Poisson variety structure {loc. tit.). 
Take a collection of loops {an} providing a pants decomposition of S. The monodromies of local 
systems over the loops a» commute under the Poisson bracket and, together with the Casimirs 
given by the monodromies around the punctures, provide an integrable system. When 5 is a 
surface with punctures it is a cluster integrable system. 

Furthermore, the Dehn twists along the loops on generate a commutative subgroup Z 39 ~ 3 in 
the modular group T(S) of the surface S, which commutes with the Hamiltonians and acts by 
cluster Poisson automorphisms. So we get a discrete cluster integrable system. 

We continue to develop these analogies in [GK]. Yet parts of this dictionary are still missing. 
In particular, it would be very interesting to construct a canonical basis - in the spirit of the 
Duality Conjectures from [FG1| - in the space of regular functions on the cluster variety X^. 



12 



1.6 Further perspectives. 

We assigned to a Newton polygon N a cluster integrable system. Here is another way to assign 
to N an integrable system: 

1. Newton polygons N, considered up to translation, parametrize smooth toric Calabi-Yao 
threefolds. Namely, the canonical bundle Kj^ of the toric surface M assigned to N is a toric 
Calabi-Yao threefold: the multiplicative group G m acts on the fibers, so the torus T x G m acts on 
Kj\f. Alternatively, it can be described as follows. Put a Newton polygon iV in the z = 1 plane 
in the lattice Z 3 , and make a cone over it. The corresponding toric threefold is a Calabi-Yao. 
Every smooth quasiprojective toric Calabi-Yao threefold is obtained this way. See |IUj for the 
connections between dimer models and toric Calabi-Yao threefolds. 

2. Any Calabi-Yao threefold Y gives rise to an JV = 2 supersymmetric Yang-Mills (N = 2 
SUYM) theory on M 4 . Namely, consider a IIB type string theory in 10-dimensional space 
Y(C) x R 4 and compactify it on Y(C). 

3. An N = 2 SUYM theory on M 4 gives rise to an algebraic integrable system - see |NSj and 
references there. 

Combining these three steps we conclude that any Newton polygon N gives rise to an al- 
gebraic integrable system. It is naturally to conjecture that it is the dimer cluster integrable 
system assigned to N. It would be interesting to prove this, and understand the role of the 
cluster structure, and in particular of quantum integrability. 

1.7 The structure of the paper. 

In Section [2] we study the set of minimal bipartite graphs T with a given Newton polygon. 

In Section [3] we define a modified partition function, and prove that it is a sum of commuting 
Hamiltonians, providing a quantum integrable system. 

In Section [J] we show that the spider move give rise to a cluster Poisson transformation. We 
show that these transformations preserve the partition functions. We also briefly recall, for the 
convenience of the reader, basic definitions of quantum cluster varieties. 

In Section [5] we study the dimer model of special kind corresponding to an arbitrary torus 
graph G, describing a resistor network on G. 

In Section [6] we construct discrete cluster integrable systems related to the dimer models on 
square grids and dimer models related to hexagonal resistor networks. 

In Section [7] we turn to the algebraic-geometric perspective of the story. As a biproduct we 
prove the independence of the Hamiltonians. 

In Section [8] we give a local calculation of the Poisson structure on the space £r, and 
generalize the story to bipartite graphs on open surfaces. 

Acknowledgments. We are very grateful to Vladimir Fock, Francois Labourie, Andy Neitzke, 
Andrey Okounkov, Dylan Thurston, Masahito Yamazaki for fruitful discussions. A.G. was 
supported by the NSF grant DMS-1059129. He is grateful to Aspen Center of Physics and IHES 
for hospitality; R.K. is supported by the NSF. 
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2 Newton polygons and minimal bipartite graphs on surfaces 

2.1 Zig-zag paths and minimal bipartite surface graphs. 

Let E = v\V2 be an oriented edge of an oriented ribbon graph. There is a notion of the left and 
right edges adjacent to E: the left (respectively the right) edge to E is the previous (respectively 
the next) for the cyclic order at V2- 

Let r be a bipartite graph on a surface S. The medial graph of T is the 4-valent graph 
with a vertex for each edge of T and an edge whenever two edges of T share a vertex and are 
consecutive in the cyclic order around that vertex. A strand is a path in the medial graph 
passing consecutively through the edges of a zig-zag path of T as in Figure EJ equivalently it is a 
path in the medial graph which goes "straight" at each (degree-4) vertex. For every edge E of T 
there are exactly two strands passing through E, see Figure El For every m-valent vertex v of T 
there is a unique m-gon formed by the parts of the strands intersecting the consecutive pairs of 
the edges incident to v, see the right picture on Figure[3l We call it the arc m-gon surrounding 
v. 




o—^ 





Figure 3: Zig-zag strands and surrounding 4-gon. 

Definition 2.1. A bipartite surface graph T is minimal if, in the universal cover T C S, strands 
have no self-intersections, and there are no parallel bigons, i.e. pairs of strands intersecting at 
two points and oriented the same way - see the left pictures on Figure [^} 





Figure 4: Minimal bipartite surface graphs have no self-intersections and parallel bigons (left); 
antiparallel bigons (right) are allowed. 



2.2 From Newton polygons to collections of admissible graphs. 

Let iV be a convex polygon in M 2 with vertices in Z 2 , considered up to a translation by a vector 
in Z 2 , called below a Newton polygon. It is described by a collection {e^} of integral primitive 
vectors with the zero sum, obtained as follows. Take the set of all integral points on the boundary 
of N. It is cyclically ordered by the counterclockwise orientation of the boundary. The vectors 
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ej are obtained by connecting each of these points with the next one. So the number of integral 
primitive vectors on the boundary coincides with the number of integral points on the boundary. 

Consider the torus T = M 2 /Z 2 . Each vector ej determines a homology class [ej] E i7i(T,Z). 
There is a unique up to translation geodesic representing this class - the projection of the vector 
ej to the torus. 

Take a family of oriented loops {«j} on the torus T in generic position such that the isotopy 
class of the loop 014 coincides with the isotopy class of the geodesic representing [ej]. The union 

G = UjOj 

is an oriented graph on the torus T. 

Definition 2.2. A graph G given by a collection of oriented loops on an oriented surface S is 
an admissible minimal graph if the following two conditions hold: 

1. Going along any of the loops on, the directions of the loops intersecting it alternate. 

2. The total number of intersection points is minimal. 

The graph G provides a decomposition of the surface S into a union of polygons Pj. The sides 
of the polygons inherit orientations from the loops. The admissibility condition 1) is equivalent 
to the following: 

V '. The sides of any of the polygons Pj are either oriented the same way, clockwise or coun- 
terclockwise, or their directions alternate. 

2.2.1 From admissible surface graphs to bipartite surface graphs. 

The admissibility condition 1') for a surface graph implies that there are three types of domains 
Pf 

1. Black domains: the sides are oriented counterclockwise. 

2. White domains: the sides are oriented clockwise. 

3. Faces: the directions of the sides alternate. 

Therefore an admissible graph G on S gives rise to a bipartite graph r on S: the black 
(white) vertices of V are the black (white) domains in S — G. Edges connecting black and white 
vertices correspond to the vertices of the graph G. Indeed, every vertex of G is 4-valent (since 
the loops OLi are in generic position), and locally looks as an intersection of two arrows. So it 
determines the black and the white domains. 

2.2.2 Minimal triple point diagrams [T] . 

A triple point diagram^ is a collection of oriented arcs in the disc with the ends at the boundary 
of the disc considered modulo isotopies, such that 

1. Three arcs cross at each intersection point. 



This is called in [TJ a connected triple point diagram. 
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2. The endpoints of the arcs are distinct points on the boundary of the disc. 

3. The orientation of the arcs induce consistent orientations on the complementary regions. 

Observe that each complimentary region to a triple point diagram is a disc. 
The orientation condition 3) is equivalent to the following: 

3'. Orientations of the arcs alternate "in" and "out" around each intersection point. 

Let n be the number of arcs. Then each of the 2n endpoints is oriented "in" or "out" , and 
these orientations alternate as we move around the boundary. Furthermore, given a collection 
of 2n alternatively oriented points on the boundary, a triple crossing diagram inducing such a 
collection provides a matching between the "in" and "out" points. 

A triple point diagram is minimal if the number of its intersection points is not bigger than 
for any other diagram inducing the same matching on the boundary. 

Theorem 2.3 (|T]). 1) In a disc with In endpoints on the boundary, all n\ matchings of "in" 
endpoints with "out" endpoints are achieved by minimal triple point diagrams. 

2) Any two minimal triple point diagrams with the same matching on the endpoints can be 
related by a sequence of 2 •<->• 2 moves, shown on Figure 





Figure 5: The 2-2 move. 



2.3 From minimal triple point diagrams to minimal bipartite surface graphs. 

Each crossing point of a triple point diagram has a canonical resolution, obtained by moving 
a little bit one of the strands passing through the point so that the emerging little triangle is 
oriented by the strands the same way as the surface - the counterclockwise resolution shown on 
Figure [6j We apply this procedure to every intersection point of a triple point diagram. 

The complimentary domains for the obtained configuration of arcs are of three possible 
types: black domains, white domains and faces, according to definitions (l)-(3) in section [2.2.21 
Namely, the black and white domains are the ones consistently oriented — counterclockwise and 
clockwise respectively — while the directions of the sides of the remaining domains, the faces, 
alternate. So we arrive at a bipartite graph associated to a minimal triple point diagram, see 
Figured Its black vertices correspond to the triple crossing points, and thus are 3-valent. The 
white vertices can have valencies two, three or higher. 

Let us shrink all edges incident to white 2-valent vertices (and remove the white vertex, 
gluing its two neighbors into a single vertex). The bipartite graph obtained can have both black 
and white vertices of valencies four or higher. 

Let us call two bipartite graphs similar if the graphs obtained by shrinking their 2-valent 
white vertices coincide. There is a canonical bijection between the faces as well as between 
the perfect matchings for similar bipartite graphs. So similar bipartite graphs have the same 
partition function, discussed in Section I3.2L 
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Figure 6: The bipartite graph obtained by counterclockwise resolution of a triple crossing dia- 
gram. 

Remark. There is another option: resolving all triple crossings clockwise. If a triple crossing 
diagram was associated to a Newton polygon N, the counterclockwise resolution is the same 
as the clockwise resolution of the diagram associated to the polygon N rotated by 180 degrees. 
Indeed, the rotation changes the orientations of all strands of the triple crossing diagram. Then 
resolving clockwise is the same as resolving counterclockwise the original triple crossing diagram. 

2.4 Boundary intervals of the Newton polygon and zig-zag paths. 

Recall that a zig-zag path has a natural orientation such that it turns right at the black vertices. 
Therefore the homology class of a zig-zag path is well defined. It is a natural invariant of a 
zig-zag path. On a minimal bipartite graph T zig-zag path has no selfintersections, and so it is 
naturally a simple loop. 

Lemma 2.4. Let N be a Newton polygon, and T a minimal bipartite graph on the torus as- 
signed to N. The homology class of a zig-zag loop is a primitive vector on a side of N, whose 
orientation is induced by the orientation of N. The cyclically ordered set of zig-zag paths with 
a given homology class is identified with the cyclically ordered set of boundary intervals on the 
corresponding side of N. So the number of zig-zag paths on T equals to the number of integral 
boundary points on N. 

Proof. It is clear from Figure [7] that zig-zag paths on the bipartite graph arising from a triple 
crossing diagram correspond to the strands of the triple crossing diagram. On the other hand 
by construction the latter correspond to the primitive vectors of the boundary of the Newton 
polygon. □ 






Figure 7: Zig-zag paths correspond to strands of the triple crossing diagram. 



2.5 Elementary transformations of bipartite surface graphs. 

Figure [8] shows an elementary transformation T± — > T2 of bipartite graphs, called a spider move. 
(There are two versions of the move depending on the color of the two interior vertices.) We will 



17 



see below that for appropriate choice of weights the spider move does not change the partition 
function for the dimer model. 





Figure 8: A spider move Ti — >■ T2 of bipartite graphs. The white vertices are attached to the 
rest of the graph; the black vertices are only attached to white vertices as shown. 

Figure [8] shows a pair elementary transformations of bipartite graphs: splitting a vertex of 
valency k > 3 and collapsing a 2-valent vertex. We will see below that it does not change the 
partition function. 

2.6 From minimal bipartite surface graphs to minimal triple point diagrams. 

A minimal bipartite graph can be transformed back into a triple point diagram in two steps. 

Step 1. By a sequence of moves inverse to the shrinking of a 2-valent white vertex, we replace 
an arbitrary bipartite graph by a bipartite graph with black vertices of valency 3, see Figure \9\ 



,0 



a. 



o 



Figure 9: Replacing a 4-valent black vertex by two 3-valent black vertices, 
splitting splits any fc-valent black vertex into k — 2 3-valent ones. 



A similar vertex 



Step 2. Take a minimal bipartite graph r on S with the black vertices of valency 3. Draw 
zig-zag strands on S assigned for all edges of F - two per edge. So each vertex of T is surrounded 
by an oriented arc m-gon. Then there is a unique way to shrink the counterclockwise triangles 
into points, getting a triple point diagram. This procedure is evidently inverse to the one we 
used to get a bipartite graph from a triple point diagram. 







Figure 10: From bipartite graphs (solid) to triple point diagrams (dashed). 
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2.7 Main result. 

Theorem 2.5. i) For any Newton polygon N there exists a minimal admissible graph G on a 
torus associated with N . It produces a minimal bipartite torus graph V associated with N . 

ii) Any two minimal bipartite graphs on a torus associated with N are related by a sequence 
of spider moves and shrinking / expanding of 2-valent vertices. 

This result is illustrated in an example in Figure HH 




Figure 11: Procedure for going from a convex integer polygon to a bipartite torus graph: 1. A 
convex integer polygon. 2. The corresponding geodesies on a torus. 3. Isotoping the geodesies to 
a triple-crossing diagram having specified order of intersection with the boundary of the square. 
4. The resulting bipartite graph. 



Proof. Let us glue a torus T from a rectangle R by identifying opposite edges. Consider a 
finite collection of oriented loops 014 on T, whose homology classes correspond to the primitive 
edge vectors of N, read counterclockwise around dN. We choose the loops so that the total 
number of intersections of any loop with an edge of R is minimal (this can be accomplished by 
starting with Euclidean geodesies, but we will isotope these around preserving the number of 
intersections with R). 

Since dN is a closed path the total homology class Yli [ai\ represented by these loops is zero. 
So the total intersection number of these loops with any edge of the rectangle R is zero. Therefore 
we can isotope the loops, by moving them one through the other, so that the intersection points 
of the loops with an edge of R alternate (and we don't add any more crossings of loops with 
dR). We can do this simultaneously for all edges of R, so that the crossings of the loops with 
OR alternate around the entire boundary of R. Then we can apply the part 1) of Theorem | 
and isotope the loops Oj into a minimal triple point diagram. 





Figure 12: Permuting the boundary intersections. 
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The components of the intersection of a loop with the interior of R will be called strands. 
Thanks to the part 2) of Theorem 12.31 an Y t wo minimal triple point diagrams whose strands 
intersect dR in the same order are related by 2 o 2 moves. We claim that by isotoping the 
loops on the torus T any two minimal triple point diagrams for the same polygon N are related 
by 2 o 2 moves. Consider the set of strands intersecting an edge E of R. Some of these strands 
may be part of the same loop, or part of isotopic loops. The relative order of intersection of 
these subsets of strands with dR is fixed. However the strands which are parts of different 
non-isotopic geodesies do not have a unique relative order along dR. We need to show that by 
isotopy and 2 -H- 2 moves we can permute the order of intersection of these sets of strands with 
E. 

Let a, c be consecutive "ingoing" strands on dR which by the alternating property necessarily 
are separated by a single "outgoing" strand b, so that a, b, c are three consecutive intersections of 
strands with dR. Suppose that the strands ending at a, c cross inside the rectangle R. We show 
that using 2 <-> 2 moves and isotopy on T we can rearrange these strands so that they intersect 
dR consecutively in the reverse order c, b, a without changing the order of the other intersections. 
In other words, we can arbitrarily permute the "ingoing" strands by doing a sequence of these 
transpositions {a, c} —> {c, a} as long as they intersect in R, and similarly for the "outgoing" 
strands. To accomplish this transposition, [T] shows that we can use 2 -(->■ 2 moves so that 
the triple point diagram inside R has the three strands at a, b, c meeting just adjacent to the 
boundary (that is, before meeting any other strands). Then one can simply isotope this triple 
crossing across the boundary which has the effect of switching the intersections order of the 
strands to c, 6, a as desired. See Figure H2] 

Now given a triple point diagram one can construct the minimal bipartite graph as above. 
Observe that there are two different 2f>2 move diagrams: the one is obtained from the other 
by reversing the strand orientations. Resolving their triple intersection points in the standard 
(counterclockwise) way, we get two elementary moves of bipartite graphs: the spider move on 
Figure [13] and the move on Figure [14] which amounts to shrinking the edge with a 2-valent 
(white) vertex, and expanding it the other way. □ 




Figure 13: This 2 <-> 2 move amounts to a spider move of bipartite graphs. 
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Figure 14: This 2o2 move amounts to collapsing the edges incident to a white 2-valent vertex 
and expanding the other way. 

3 Integrability of the dimer model 

3.1 Edge weights and line bundles with connections on a graph 

A vector bundle with connection on a graph T is the following data: we assign to each vertex v of 
r a vector space V v , and to each oriented edge E from a vertex v% to a vertex V2 an isomorphism 

IE ■ V Vl -s- V V2 

called parallel transport, so that i- E = i E , where — E is the edge E with the opposite orienta- 
tion. The monodromy of the connection along a loop is the composition of the parallel transports 
around the loop. Two vector bundle with connections {V,%e}-, {V,i' E } are isomorphic (gauge 
equivalent) if there are isomorphisms ip v : V v — > Vy such that i' E = tp V2 o i E o -0" 1 for all edges 

E = V\V2- 

The moduli space £r is the space of line bundles with connections on a graph T modulo 
isomorphisms. Here is another description of the moduli space C-p. Choose a non-zero vector e v 
in each space V v . Then a line bundle with connection on T amounts to a collection of non-zero 
numbers, the edge weights a E , assigned to the oriented edges E of the graph, so that a- E = a E : 

i E e Vl =a E e V2 , E = v\ -+ v 2 . 

Another choice of base vectors e' v = b v e v gives rise to a different set of edge weights a' E , related 
to the original one by a gauge transformation: a' E = b v a E b~, . So 

edge weights 



gauge transformations 

Given a trivializations {e„} of a line bundle with connection {V,i E }, we can encode the 
connection by the Kasteleyn operator K, see Section [1.4.21 which in this case is just a matrix. 
However the determinant det(K) is not gauge invariant, and so is not a function on the moduli 
space Cr- 

3.2 The partition function and the Hamiltonians 

3.2.1 The weight of a dimer cover 

A dimer cover of a bipartite graph T, also known as a perfect matching, is a collection of edges 
of r such each vertex is the endpoint of a unique edge in the dimer cover. 
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Given a connection, a dimer cover M on T provides an L-valued function W(M) on £r> 
called the weight of M. For an edge E = bw, let %e be the parallel transport along E of L. Then 
iE € V£ <S)V W . The product of these vectors over the edges of M is the weight: 

W{M) €0(£ r )®L. 

Precisely, choose an order of the edges E±,..., E^ of the matching. Write %e = a* E <g> frgj. Then 
VF(M) := a^ A ... f\a* E ® bE x A • • • A bE N - It does not change when we alter the order of edges. 
The following result, due to Kasteleyn (in different language), is well known. 

Theorem 3.1 ([Kast]). We have 

detK = ^2sgn(M)W{M) eO(£ r )»L (17) 

M 

where the sum is over dimer covers and the sign sgn(M) only depends on the choice of k and 
the homology class modulo 2 of M (see definition below). 

Definition 3.2. The partition function is Vv = det-K". 

The partition function depends on our choice of k, which we fixed once and for all. 

3.2.2 The homology class of a dimer cover 

A dimer cover M on T provides a function 

com : {Edges of T} ^ {0,1}. 
Its value on an edge is 1 if and only if the edge belongs to M . We view it as a 1-chain 

[M] := -J2"m(E)[E], 

E 

where the sum is over all edges, and [E] is the edge E oriented from black to white. The 
condition that M is a dimer cover means that 



d[M] = Y,sm(v)[v] 



where the sum is over all vertices v of I, and sgn(u) = 1 for the black vertex v, and sgn(u) = — 1 
for the white. 

If M\,M2 are dimer covers, [Mi] — [M2] is a 1-cycle, and thus each dimer cover defines an 
"affine" homology class. It is more convenient to define an absolute homology class for each 
dimer cover by fixing a chain $ such that [M] — $ is a cycle (for all M). We do this in Section 
I3T2~31 

Remark. A dimer cover M2 trivialises the determinant line L by providing isomorphisms 
Vb — ► Vw for each edge E = bw. In Section F3.2.3I we consider "fractional" dimer covers provided 
by weighted zig-zag paths, and use them a similar way to trivialize L. 
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3.2.3 The a-deformed partition function 

For every edge of an oriented ribbon graph there are exactly two zig-zag paths containing the 
edge (one in which the edge is traversed in either orientation). Let E be an edge of V. Denote 
by z r (respectively z) the zig-zag path containing E in which the white vertex of E precedes 
the black (respectively follows the black). Let us consider a map 

a : {zig-zag paths on T} — > R/Z. (18) 

It provides a function 

tp a = <p : {Edges of T} — > R, <p(E) := a r - a L . (19) 

where a r (respectively a{) is the evaluation of the map a on the zig-zag path z r (respectively 
z) crossing the edge E, see Figure [T5| and a r — a,\ denotes the length of the counterclockwise 
arc on S l = R/Z from ai to a r . 



Figure 15: The weight is <p(E) = a r — a.\. 

The strands, and thus the oriented zig-zag paths, are in bijection with the set {ei, . . . , e n } 
of primitive edges of the boundary dN of N, and so inherit a circular order from ON. 

Theorem 3.3. Let X be the set of circular- order-preserving maps 

{ei,...,e n }^R/Z. 
For any map a € X, the function ip a satisfies 

dp a = £)sgn(«)[t;]. (20) 

V 

The result is more generally true for arbitrary triple crossing diagrams in a disk. 

Proof. We need to show that dip a = 1, respectively —1, at a black, respectively white, vertex. For 
this it suffices to show that at each vertex, the set of adjacent strands (the strands intersecting 
edges emanating from the vertex) has the same circular order as the circular order of the edges 
containing the vertex. This follows from Lemma 13.41 below. □ 

Lemma 3.4. Given a triple- crossing diagram in a disk, at any triple- crossing the three inter- 
secting strands occur in the same circular order as their ingoing endpoints on the boundary of 
the disk. 
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Proof. The proof is by induction on the number of strands. Order the strands according to the 
counterclockwise circular order of their ingoing endpoints. We recall the construction of a triple 
crossing diagram of [T]: take a strand with the property that on one side, no other strand has 
both its endpoints. One can draw this strand following the boundary as illustrated in Figure [TBI 
crossing the intermediate strands in pairs. Each of these crossings is positive. The remaining 
region again has the property that the strands alternate orientation along the boundary and so 
can be completed to a triple crossing diagram with one fewer strands by induction. Moreover 
the ingoing strands of the new region retain their original circular order. □ 




Figure 16: 

By taking limits in which all points a(e^) collide (but preserve their circular order) there are 
circular-order-preserving maps a for which ip a takes integer values <p a € {0, 1}. 
Given a map a provided by Theorem 13.31 there is a 1-chain 

$ « = Yj ^ a ( E ) [^ with d ® a = J2 sgn ( y ) M • 

E v 

It is an integral 1-chain if cj) a is integer-valued. Thus [M] — & a is an integral cycle. It provides 
a monomial function on £p 

GJ a (M):£ T — ^C*, (21) 

given by the monodromy around the cycle [M] — & a . It is the complex character of the group 
Ct corresponding to the homology class [M] — $ Q via the isomorphism (pQ). 

Definition 3.5. Let T be a bipartite surface graph. Let a be a map satisfying A20\). The a- 
deformed partition function V a of the dimer model on V is a function on the space Cy given as 



the sum of monomials Lo a {M) over all dimer covers M ofT, with the same signs as in J77[ ): 



P Q :=^sgn(MK(M). (22) 

M 

3.2.4 The Hamiltonians 

Let [M] a G Hx(T, Z) be the homology class of the cycle [M] — $> a . Recall the projection 
it : i?i(r,Z) — > Hi(S,Z). Given a homology class a € Hi(S,7<), take the part of the sum 
over the matchings M such that the projection of [M] a to Hi(S, Z) is the class a: 

H a - a := Y sgn(M)u a {M). 

M-.Tr([M] a )=a 
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By Theorem 13. II the sign sgn(M) depends only on a, and the partition function is a sum 

V a = Y, sgn(a)tf Q;a . (23) 

aeHi(S,Z) 



Definition 3.6. The functions H a - a in /i23\) are the Hamiltonians of the dimer system. 

Remark. A different map a' leads to another collection of Hamiltonians {H a /. a } which differ 
from {H a;a } by a common monomial factor, which lies in the center of the Poisson algebra. 
Indeed, the monodromy along any zig-zag loop is a Casimir element. Therefore the Hamiltonian 
flows do not depend on the choice of a, and on each symplectic leaf any two collections of 
Hamiltonians differ by a common constant. 

A coordinate expression. Let zx,...,Z2 g be loops on T whose homology classes generate 
Hi(S,Z). Let ^nj[zj] be the homology class of the cycle [M] — $ a . Then there is a function 
ipM '■ {Faces of T} — > M, which we interpret as a 2-chain, such that 

dlp M = [M}-$ a -^2 n A z iY 
So one has 



V a :=Y±U W F } Il^- 



Am{F) 

vv 
M F i 

Each Hamiltonian is, up to a sign, a sum of monomials with coefficients +1 in variables Wf, Zi 
with a given value of (n\, ...,ri2g)- If S is the torus, P a is a polynomial in z%, %2- 

3.3 Complete integrability 

A minimal bipartite graph on a torus provides an algebraic quantum integrable system: 

Theorem 3.7. Let V be a minimal bipartite graph on a torus T. Then 
i) The Hamiltonians H a - a commute under the Poisson bracket on £r- 
ii) The Hamiltonians are independent and their number is the half of the dimension of the 

generic symplectic leaf. 

Proof, i) Take a pair of matchings (Mi,M2) on V. Let us assign to them another pair of 
matchings (M\,M2) on T. Observe that [Mi] — [M2] is a 1-cycle. It is a disjoint union of 

1. homologically trivial loops, 

2. homologically non-trivial loops, 

3. edges shared by both matchings. 

For every edge E of each homologically trivial loop we switch the label of E: if E belongs to a 
matching M\ (respectively M2), we declare that it will belong to a matching M2 (respectively 
Mi). For all other edges we keep their labels intact. Denote by \Xi and Jli the homology classes 

of the cycles [Mi] - $ a and [M] - $ a . Clearly (M, M 2 ) = (M 1 ,M 2 ). We use a shorthand e(*, *) 
for £(*,*) § r . 
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Lemma 3.8. One has 

e(Mi>M2) + e(Mi,M2) = 0. (24) 

Proof. We use a local calculation of the Poisson bracket on £-p given in Section Let us write 
(j24[) as a sum of the local contributions corresponding to the vertices u of T, and break it in 
three pieces as follows: 

£(AJi,M2) + e(Mi,M2) = (Yl + ^2 + ^)(^u(M1>M2) + 5u(M1,M2)J- 

i>e£i ?>e£2 ^££3 

Here £i (respectively £2 and £"3) is the set of vertices which belongs to the homologically trivial 
loops (respectively homologically non-trivial loops, double edges). 

Then for any v G £ 1 we have 6 v (fj,x, ^2) + ^(^1,^2) = since the elements l v (fj,x) € A^ 
provided by //, by the very construction coincides with the element l v (Ji2), and similarly l v (/j,2) = 
l v {Jli). So the first sum is zero since it is a sum of zeros. 

The third sum is also a sum of zeros. Indeed, the contribution of the matchings M\ and M2 
at every vertex v € £ 3 is the same, so the local pairing vanishes by the skew symmetry. 

It remains to prove that 

J~] ( S v (fii, 1x2) + S v (Jli, Jj 2 )) =0. 
ve£ 2 

Let 7 be an oriented loop on T. We define the bending 6^(7; </?) of the function cp at a vertex 
v of 7 as 

6,(7;^)= 2 v?(£)- ^ ^)el 

Here -R„ (respectively L^) is the set of all edges sharing the vertex v which are on the right 
(respectively left) of the path 7, following the direction of the path 7. Lemma 13.81 now follows 
from Lemma 13.91 below. □ 



Recall the set X of circular-order-preserving maps from the zig-zag paths to K/Z. 

Lemma 3.9. There is an a € X such that the corresponding function ip a satisfies, for any 
simple topologically nontrivial loop 7, 

Y^bvir^a) = 0. 

ve-y 

Proof. On a Euclidean torus, a simple closed polygonal curve which is not null-homotopic has 
the property that its net curvature is zero, where the net curvature is the total amount of left 
turning minus right turning. This fact can be seen by taking an isotopy of the curve to a straight 
geodesic through a family of simple polygonal curves; the total curvature is an integer multiple 
of 2-7T but changes continuously under the isotopy, and for a geodesic it is zero. 

A map a € X allows us to draw our graph T on a Euclidean torus in such a way that 
6^(7, (p a ) is (with a small proviso) the net curvature of 7 at vertex v. This is accomplished as 
follows. To each edge E, recall that a r , a\ are the angles of the strands crossing E. We define a 
I-cochain / by f{E) = e l ° l + e tar . We have df = 0, since strands cross two edges of every face, 
once with each orientation. Thus there is a multivalued function F on vertices of T, which is 
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locally well defined but has additive periods u)\,U2 for the generators of H\{S). Extending F 
linearly over the edges gives a map from T to a Euclidean torus C/(wi, 0*2). Locally, to each edge 
we are associating a Euclidean rhombus with vertices z, z + e tOLl , z + e iai + e iar , z + e tar , where z 
is determined by F; the edge -E 1 runs from z to z + e* Q ' +e JQr , and the vertices z + e iai , z + e l ° r are 
the centers of the faces adjacent to E. These rhombi are glued edge to edge to form a Euclidean 
torus (possibly with conical singularities). 

If for each E the angle (p(E) = a r — ai is less than ir, then the rhombi are positively oriented 
and, by Lemma 13.31 the total angle at each vertex is 2tt and our Euclidean surface is smooth, 
that is, has no singularities. If a rhombus for edge E, connecting vertices v,v', has angle larger 
than 7T, by Lemma 13.31 it will be the only rhombus at v or v' with this property. Hence the 
singularities in our Euclidean structure are folds, as shown in Figure [T7"l and are separated from 
each other. 

It remains to check that if 7 contains a vertex v which has a singularity then it is still true 
that 6^(7, ip) measures the net curvature of 7 at v. If 7 passes through v but not along edge 
E this is trivially true. In the remaining case, 7 contains edge E = vv': suppose 7 contains 
consecutive vertices vi,v, v' , V2 in order. Then it is straightforward to check that the sum of the 
curvature contributions at v and v' equals the change in Euclidean angle of the edge v\v to the 
edge v'v2- In all cases therefore b v is the net curvature of 7, which is zero. □ 




Figure 17: 

End of proof of Theorem \3. 1$ ). The part i) of Theorem 13.71 follows immediately from Lemma 
3.81 Indeed, the sum calculating the Poisson bracket of two Hamiltonians can be presented as a 
sum of contributions (1241). □ 



The part ii) of Theorem 13.71 is proved in Section 13.61 below. 

3.4 The matching polytope 

Let f2 be the matching polytope of T. By definition this is the subset 

ft C [0,l]i ed g csofr l 

consisting of functions uj on edges of V with values in [0, 1] such that the sum of the values at 
every vertex is 1. The subset ft is a convex polyhedron. 

Lemma 3.10 ( |LPj ). The vertices of ft are the perfect matchings ofT. 
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Proof. To make the paper self-contained, we give a proof. Let w be a vertex of f2. Let us show 
that for every vertex v of T there is just one edge E incident to v with uj(E) ^ 0. Indeed, 
otherwise there is a vertex of v\ incident to edges E\,Ei with the weights in (0, 1). Let v-i be 
another vertex of E^. It is incident to another edge E3 with the weight in (0, 1) - otherwise the 
sum of the weights at V2 will be less then 1. Let V3 be another edge of E3, and so on. We get a 
closed path on V whose edges have the weights in (0, 1). Let us add e to the weight of each of 
these edges, where \e\ is very small, and the signs alternate. The sign condition does not lead 
to a problem since T is a bipartite graph. Then we still get an element u e € 0. So uj is inside 
of a segment in Q,. Thus it is not a vertex of fi. 

It follows that the values of uj at the edges are either 1 or 0. The edges where the value is 1 
give a matching. □ 

Elements of £1 are called fractional matchings. 

Lemma 3.11. Let Y be a minimal bipartite graph on a torus. Then Y has a perfect matching. 

Proof. It suffices to show that f2 is nonempty. Indeed, by Lemma 13.101 any vertex of the poly- 
hedron fi provides a perfect matching. Any minimal bipartite graph Y is obtained from some 
polygon N. The function ip a from Lemma 13.31 is then a fractional matching by (|20p . □ 

Let ujq € f2 be a fixed fractional matching. Then for any u G Q the difference [uj] — [ujq] is 
a 1-cycle on T. Taking its homology class we get a map 0, — >■ Hi(T,M). Its image is a convex 
polyhedron [O] C Hi(T,M). Its projection under the canonical map ir : iTi(r,R) — > Hi(T,M) is 
a convex polygon 

ir[n] c Hi(T,R) ^R 2 , 

which depends on ujq only through a translation. 

Recall the space X of order-preserving maps from {e\, ... ,e n } to R/Z. The image of X 
under the map a —> <p(a) is a subpolyhedron Qq C 0. It is in fact a simplex: the simplex of 
order-preserving maps from the strands to R/Z, modulo rotation. 

Let on £ X be the map which sends all E{ to and assigns 1 to the arc between E% and -Bj+i 
and to all other arcs. The matching (p(cti) maps to a vertex of Qq, and every vertex of ^o is 
of this type. 

3.5 The Newton polygon 

Theorem 3.12. Let T be a minimal bipartite graph with polygon N. Then the Newton polygon 
of the partition function P a {z\, z%) for T is N (up to translation). 

Proof. Set wo = 3> a . The Newton polygon of P a (z\,Z2) is a translate of 7r[f2]. Indeed, by 
definition the coefficient in P of z 1 ^ z % <£ is the sum of weights of dimer covers M of Y for which 
the homology class [M] — $ Q satisfies [M] — $ Q = ii[zx] + ^2^2]- 

We can define the intersection pairing of a chain with a zig-zag path Z as follows. If C = 
^eCe[E] is a chain, the intersection pairing of C with Z is (C, [Z]) = Y^egz^Ce where the 
sign is alternating, and depends to whether the orientation of the edge of Z agrees with the 
white-to-black orientation or disagrees with it. In other words, this is the intersection pairing 
of C with the strand corresponding to Z. 
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We define another polygon Nr C i?i(T,R) by a set of inequalities 

N r := J 7 £ #i(T, E) | (7, [Z]> < ^ " ([wo], [£]> for every zig-zag path z| . (25) 

We claim that N, N^ and ir[Q] are translates of each other. 

We first show that ir[Q] C Nr- Let uj £ CI. The sum of u(E) over the edges Ein a zig-zag 
path Z is at most equal to l(Z)/2, since this is the number of white vertices adjacent to Z (each 
white vertex adjacent to Z has total flow at most 1 across Z). Thus 

<[ W - «*,], [z]) = <[«], [z]> - (M, [z]> < ^ - (N], [z]> 

which implies that 7r[fi] C iVr. 

We now show that the edges of N and iVr coincide. Take 7$ = [ui(ai) — ujq\. This satisfies 
the inequalities (f25j) and two of these are equalities: those for Z% and Zj+i, because in «j every 
white vertex adjacent to Zj is matched to a black vertex on the other side of Zj, and similarly 
for Zj + i. Note moreover that 7j + i = 7, + [Zj] so we see that [Zj] = e, is an edge of iVr- Thus 
Nr has sides consisting of the ej in order, which proves that N and Nr are translates. Note 
moreover that each corner of iVr is in the image of ir[Q] (in fact in the image of 7r[f2o])- Thus 
tt[0] =N r . ' " D 

3.6 Counting triple crossings and faces 

From Thurston (Tj, the number of triple crossings in any triple crossing diagram in a disk is a 
function solely of the strand order along the boundary, and can be obtained as follows. Take 
a linear order of the strand endpoints which is consistent with the circular order. Then we say 
two strands ac and bd, directed from a to c and b to d respectively, have a parallel crossing if 
their endpoints are in the order a < b < c < d or d < c < b < a. 

Lemma 3.13 ([T]). In a minimal triple- crossing diagram in a disk, the number of triple cross- 
ings is the number of pairs of strands with a parallel crossing. 

Now take a triple crossing diagram on the unit square [0, l] 2 obtained from a convex integral 
polygon N as in Section [2j Let ei, . . . ,e n £ Z 2 be the primitive sides of N. 

Lemma 3.14. The number of triple crossings is twice the area of N. 

Proof. A minimal triple crossing diagram can be obtained by starting with a configuration of 
straight geodesies on the torus [0, l] 2 , one for each e*, and then isotoping them into position. 

We count the number of parallel crossings in a standard configuration of the geodesies ej and 
then notice that this number is invariant when we isotope the strands to their final position. 

Let r £ R 2 be a vector with positive coordinates having irrational ratio, and let v m i n ,v max 
be respectively the unique vertices of N which minimize, respectively maximize v ■ t. Let Nr 
be the set of sides of N between v m i n and v max in eclw order and Nl be the set of sides of N 
between v max and v m i n in eclw order. Then after a circular re-indexing, say that e%, . . . , e& £ Nr 
and e/c+i,. . . ,e„ £ N L . 

Let p = et £ [0, l) 2 for some sufficiently small e. Start by drawing all geodesies ej in the 
torus through p. No geodesic contains the origin by irrationality of r; order the boundary points 
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on the square cclw starting from the origin. Then ej and ej have a parallel crossing at p if and 
only if ej and ej are either both in Nr or both in Nl. Moreover if p is sufficiently close to the 
origin all other crossings of ej and ej will be parallel or non-parallel according to the crossing at 
p. 

Thus the number of parallel crossing strands in total is 

y] |ejAej[+ y |ejAej|. 

However the area of the part of TV below the line through v m i n and v max is 1/2 of 

ei Ae 2 + (ei + e 2 ) A e 3 H h (ei H he fc _i)AeA;= ^ |ejAej|, 

i<j<j<fc 

as can be seen by triangulating this polygon using extra edges from v m i n to each other v. 
Similarly for the upper part of N. 

Now as we isotope the geodesies ej, without letting a strand cross the origin, it is easy to see 
that the number of parallel crossings does not change. 

□ 

Let us denote by B and W the number of black and white faces for a triple point diagram 
graph G on the torus arising from a Newton polygon N . Let V be the number of vertices of G, 
that is the number of the triple crossings. 

Proposition 3.15. One has 

B = W = V = number of the triple crossings = 2 (area of N). 

Proof. Euler's formula for the triple point diagram G gives V — E + F = 0. Counting flags 
(vertex, edge) we get QV = IE. Therefore F = 2V . 

Lemma 13 . 1 1 1 implies that B = W . Indeed, the white (respectively black) domains of a triple 
point diagram G became the white (respectively black) vertices of the bipartite graph I\ Thus 
F = B + W implies B = W = V. Lemma 13.141 implies the last statement. 

□ 

End of the proof of Theorem 13.71 Let i(N) and e(N) be the number of interior and 
exterior points for the Newton polygon N. Pick's formula plus Proposition 13.151 gives 

2i(N) + e(N) -2 = 2(area of N) = B. (26) 

Observe that i(N) is the number of Hamiltonians, and e(N) is the number of Casimirs plus one, 
that is the dimension of center of the Poisson algebra O(Cr) minus one. On the other hand, B 
is the number of faces minus one, and ikHi(S, Z) = 2. So B + 1 = dinnCr- So (|26p implies that 



2(number of Hamiltonians) = dim£r — dim I Center O(Ct) 

The independence of the Hamiltonians is established in Corollary |7.5l below. which is independent 
of the proof of the theorem. 

iii) Quantum integrability. The claim that the g-deformed Hamiltonians commute follows 
immediately from the combinatorial proof given above. This completes the proof of Theorem 
13771 □ 
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4 Dimers, cluster Poisson transformations, and quantum inte- 
grability 

4.1 Dimer models on surfaces and cluster Poisson varieties. 

4.1.1 Gluing the conjugated surfaces under spider moves 

Given a bipartite ribbon graph T, let Sr be the corresponding black conjugated surface. Since 
the conjugated graph Tj, coincides with T as a graph, there are canonical isomorphisms of lattices: 



H 1 (S r ,Z) = H 1 (r,Z) = H 1 (T > Z). 



(27) 



We denote by Ar the lattice (|27p . It is equipped with the skew symmetric integral bilinear form 
(*, *) given by the intersection pairing on Sr- The boundaries dF of the faces F of T, whose 
orientations are induced by the orientation of S, give rise to a collection of cycles {7_f} on Sr- 
These cycles plus the generators ai, ..., a<i g of Hi(S, Z) generate Ar, and satisfy the only relation 
Ef7f = 0. 

A spider move is a transformation of a bipartite graph shown in dashed lines on Figure [T8l 
We show on the same figure the transformation of the accompanying zig-zag paths (solid) , and 
the unique cyclic order (1,3,2,4) of the zig-zag paths compatible with the cyclic orders of the 
triples of the zig-zag paths at the black vertices. These cyclic orders are the same for both 
graphs. 




I i 



Figure 18: A spider move and zig zag paths. 

A spider move of bipartite ribbon graphs s : T — > V is determined by a quadrilateral face 
F of r. Denote by Fi , F2 , -F3 , F4 the four faces of T sharing edges with F. Denote by Sr and 
Sr' the corresponding conjugated surfaces. The spider move does not affect the graph V outside 
of a little domain U containing the face F and the two extra edges emanating from the two 
black vertices of F. So we identify the conjugated surfaces Sr and Spi outside of U. We want 
to extend this identification to a homeomorphism of surfaces 



's : Sr — > St> ■ 



(28) 



Lemma 4.1. There is a unique up to isotopy homeomorphism &28\) identifying the boundary 
loops. Let j' G := s"(7g). Then one has 



l' F 



-IF, 



7Fi = 7Fi + IF, lF 2 = 7F 2 , 7f 3 = 7F 3 + IF, 7f 4 = 7F 4 - (29) 



and other classes jg remain intact. 
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Figure 19: Identifying the conjugated surfaces according to a spider move. 

Proof. Figures [19] and [20] show the ribbons of the conjugated graphs T and V inside of the 
domain U; in Figure [20] the ribbon graph is shown in its untwisted view. These domains are 
annuli. 

Pairs of parallel zig-zag paths together with a pair of little segments connecting them form 
the boundary of the annuli. 

The "before" and "after" annuli have the same boundary and so can be identified modulo a 
half of a Dehn twist, which we perform as in Figure [20l So we get the first claim of the Lemma. 

Looking at the loop jf (red dotted loop on Figure I20J) and its image after the spider move 
we see that 7^ = —7^. Similarly the remaining formulas of (|29p follow from considering the 
corresponding zig-zag paths on Figure [20] □ 

4.1.2 Seeds, mutations, conjugated surfaces and spider moves 
Definition 4.2. A seed s is a datum ( A, {e{\, (*, *) I , where 

• A is a lattice, i.e. a free abelian group, {ei} is a collection of non-zero vectors A; 

• (*, *) a skew- symmetric %-valued bilinear form on A. 

Definition 4.3. A mutation of a seed s in the direction of a basis vector e k is a new seed s. 
The lattice and the form of the seed s are the same as of s. The collection of vectors {ej} of s 
is defined by applying a half reflection at e k map to the old one: 



e,: := 



ei + (ei,e k ) + e k if i^k 
-e k if i = k. 



a + := max(a, 0). 



(30) 



Definitions 14.21 and 14.31 generalize slightly the standard definitions of seeds and seed muta- 
tions, see Sections 1.2.1-1.2.4 of [FG1]. Unlike loc. cit, we do not demand neither that the 
vectors e, are linearly independent, nor that they span the lattice Ao 

Here is our main example. 



'We simplify the discussion by considering only the "simply-laced" case, with the "multipliers" di 
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Figure 20: Identifying the conjugated surfaces according to a spider move, untwisted view. The 
red dotted curves are the curves jp and j' F . The blue dotted curves are part of jp 1 and j' F . 

Definition 4.4. The seed sr assigned to a bipartite graph T on a genus g surface S is 

s r := (A r ,(*,*)p,{7F})- 

Lemma 4.5. The spider move T — > V centered at a face F leads to a mutation of seeds 
sr — > Sr*' in the direction of the vector jp . 

Proof. The only non-zero values (7f,7g)p are 

{lF,lF 1 )f = l, (7F,7F 2 )f = -1, (7F,7F 3 )f = 1 > (7F,7F 4 )f = -l (31) 

Thus formalae (129p from Lemma 14. II are just equivalent to the mutation formulae (130|) . D 

Therefore the geometric objects - the conjugated surfaces and the spider moves - are encoded, 
respectively, by the corresponding seeds and their mutations. 

4.1.3 Quantum cluster transformations [FG1] 

A lattice A with a skew-symmetric Z-valued bilinear form (*,*) gives rise to a quantum torus 
*-algebra Ta- The algebra Ta has a basis {X v } over the ring Z,[q,q~ 1 ] parametrized by the 
vectors v of the lattice A. The multiplication is given by 

n -(vi,V2) y Y — Y 

There is an involutive antiautomorphism 

*:T A ^T A , *(X v )=X v , *(q)=q-\ 

Let us choose an order of the basis {vi} of A. Then 



n n 



Xv = g -E«,«i«yKi*) Y[X? , v = ^OiVi, (32) 



t=i »=i 
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Consider the following formal power series, called the q- exponential: 

It is characterized up to a constant by the difference relation 

^ q {q 2 x) = {l + qx)^ q {x). (34) 

Denote by Ta the non-commutative fraction field of the algebra Ta. The quantum mutation 
map n% k in the direction of a vector e^ is an automorphism of the skew field Ta given by the 
conjugation by & q (X ek ): 

4 := Ad* q{Xek) : T A — > T A . (35) 

Although *& q (X ek ) is not a rational function, this is a rational map. 

4.1.4 Cluster Poisson varieties [FGlj 

Cluster Poisson varieties, also known as cluster ^-varieties [FGlj . are geometric objects which 
are in duality to cluster algebras of Fomin-Zelevinsky [FZIJ. 

We assign to a seed s a complex algebraic torus X\ := Hom(A, C*), called the seed torus. A 
vector v E A gives rise to a function X v on the torus obtained by evaluation at v. A mutation 
in the direction of a vector e^ gives rise to a birational isomorphism [i ek : X\ — > X\, obtained 
by specialising q = 1 in (|35[) . It acts on the function X v by 

£ h :X v *-+X v (l + X eh )-te> e >\ (36) 

The Poisson bracket is the quasiclassical limit of the commutator: 

{X n ,X V2 } := lim[X !)1 ,X y2 ]/2(g< - 1) = (v 1 ,v 2 )X Vl X V2 . 

q-+l 

Thus (|36p is a Poisson map. An isomorphism of seeds leads to a Poisson isomorphism of the seed 
tori. Compositions of seed mutations and seed isomorphisms are seed cluster transformations. 
So a seed cluster transformation gives rise to a Poisson birational isomorphism of the seed tori. 
A seed cluster transformation s — y s acting trivially on the functions X v is a trivial seed cluster 
transformation. Consider all seeds obtained from a given seed s by a sequence of mutations. 
We assign to these seeds their seed torus, and glue pairs of tori related by mutations into a 
space X by using birational isomorphisms (|36p . Trivial seed cluster transformations give rise 
to automorphisms of X. Taking the quotient of X by these automorphisms, we get a cluster 
Poisson variety X[j 

Non-trivial cluster seed transformations s — > s lead to automorphisms of the X, which form 
the cluster modular group. We will see such examples in Section 6. 

A cluster Poisson variety is nothing else but a collection of seed tori X\ glued via the cluster 
transformations. The algebra of regular functions 0{X\) on a cluster torus is an algebra of 
Laurent polynomials. So a cluster Poisson variety is described by a collection of algebras of 
Laurent polynomials and birational automorphisms (j36[) providing the gluing transformations. 



3 It could be a non-separated space rather then a variety. 
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The algebra O(X) of regular functions on a cluster Poisson variety X can be identified with the 
subalgebra of the algebra O(X^) consisting of the elements which remains Laurent polynomials 
after any cluster transformation. 

This suggests that one should think about the quantum torus algebra Ta as of the algebra 
of regular functions on a non-commutative quantum torus, and that these algebras together 
with the mutation maps (|35p should be thought of as a quantum cluster variety. We define 
the algebra O q {X) of regular functions on a quantum cluster variety as the subalgebra of the 
quantum torus algebra Ta consisting of all elements which remain Laurent polynomials after 
any quantum cluster transformation. 

4.1.5 Gluing the moduli spaces Cr according to a spider move 

Instead of using the set of all seeds related to an initial seed s by seed cluster transformations, 
we can pick any subset of this set, and define a Poisson space by gluing the corresponding cluster 
tori. The cluster Poisson spaces which we need in our paper are defined this way: we use only 
finitely many seeds corresponding to minimal bipartite graphs associated with a given Newton 
polygon. 

4.1.6 Traditional cluster Poisson variety related to the dimer model 

Traditional cluster Poisson varieties are the ones defined by using seeds where the vectors {ej} 
form a basis of the lattice. Let 

T s := Ht&C*) = (C*) 2g . 

It is the group of line bundles with flat connections on S. The embedding T ^ S provides a 
free Poisson action of the group T$ on the moduli space £r- 

The monodromy of a line bundle on T around a face F is a function Wp on C-p- The product 
of all Wp's is 1. Define 

X' = ('(T*Uf aces °f ^} 

Let W : <Yp — > C* be the product of all coordinates. Assigning to a line bundle on T its 
monodromies around the faces we get a canonical isomorphism 



(£ r /T s )^Ker(x^ 



The birational Poisson transformations ([7]) are Tg-equivariant. So to describe them on the 
quotient is just the same as to describe their action on the face weights. 

Lemma 4.6. The algebraic tori {X^} parametrized by equivalent minimal bipartite graphs T on 
S are glued by cluster Poisson transformations {7|) into a cluster Poisson variety X' N . The maps 
([?[) preserve the function W. So there is a hypersurface X^ C X' N given by the equation W = 1. 
The group T$ acts freely on Xn . There is a Poisson isomorphism 

X N /T S ^ X%. (37) 

Proof. We just need to show that the product of the face weights is preserved by the spider 
move cluster trabnsformations. This is clear from formulas (1381) below. □ 
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4.2 Spider moves and partition functions 



We have seen that elementary transformations r — >• V give rise to rational maps between the 
moduli spaces £r —* £>v- m this Section we show that they do not change the partition functions 
for the dimer models related to the graphs V and V . Moreover, we show that the rational maps 
are characterised by the condition that they preserve the partition functions. 

Theorem 4.7. Given a spider move, there is a unique rational transformation of the face 
weights preserving the modified partition function V a . This transformation is a cluster Poisson 
transformation, see Figure [2l[ 



(38) 
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Figure 21: A spider move and a mutation of the face weights. 



Proof. Let us show first that formula (|38|) is a cluster Poisson transformation. Precisely: 



Lemma 4.8. Formula \38\) describes the cluster Poisson mutation centered at jf, where F is 
the center of the spider move. 



Proof. Let us compare formula (|38p with the formula for the cluster transformation at 7_p. For 
this purpose, let us write the action of the mutation at the vector et on the functions X e > in 

i 

terms of the original functions X &i : 



/V : x e 



X 



-i 



if i = k, 



\ Xi(l + x- B8n(ei ' efc) )-( e *' e *> if i^k. 
Evidently, thanks to ()3ip , formula (j38j) is a specilaization of the general formula ([39 



(39) 
D 



Proof of Theorem \4-7\ We use the edge weights, where the notation are as on Figure 
The equality of the modified partition functions reads as follows: 

Pr_ P^ 



(40) 



where Pp is the partition function of the dimer model on the bipartite graph T defined using 
the edge weights, and 

M a := Y[ weight(£)^. 

edges E 
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Here (p(E) = a r — a\ is determined by the function a, see (USD, CEHD, an d Figure [T31 Observe 
that the maps a for the graphs V and r" are the same maps. 

The face weights are calculated via the edge weights as follows, where W* are factors common 
for both graphs T and r", and the edge weights are as on Figure [ 



dc b ci & f 

W = —, W 1 = -Wf, W 2 = -W 2 * ,W 3 = -W£, W A = J -Wl. 
be a c j a 



W 



AC 
BD' 



W[ 



A 



Wf, Wn 



1 - 



B 
E 



w* 2 ,w£ = §wz, wi = j\vi 




Figure 22: A spider move and a mutation of the edge weights. 

Let us number the external vertices on the graphs shown on Figure [22] by 1,2,3,4. We 
present the graphs as unions of the inside and outside graphs, where the inside graphs are the 
ones shown on Fig[22j Given a pair of vertices i,j S {1,2,3,4}, consider all perfect matchings 
where the vertices i and j are matched to the outside vertices. Their contribution to the partition 
function is the product of the total weight of the perfect matchings for the outside graph by the 
weight pij of the perfect matchings for the inside graph. Similarly p^ ■ for V. One should have 



Pij — Cpij, 



(41) 



with the coordinates (pi 2 ■ P13 ■ ■■■ ■ P34) and (p' l2 '■ P13 '■ ■■■ ■ p'34) coincide. 



i.e. the points in 
One has 

P12 = bf, pi 3 = be + dc, pu = cf, p 2 3 = ad, p 2 4 = af, p 3 4 



ae. 



P12 



CF, p' 13 = EF, p' u = DE, p' 23 = BF, p' M = AC + BD, 



Pu 



AE. 



(42) 
(43) 



For example, there are exactly two matchings of the inside graph which do not use the vertices 
1, 3, and their total weight is ^13 = be + dc. 

Then for every pair (ij) and {kl) there is a condition 



Pii 
Pkl 



Al 



Let us investigate its consequences. We deduce that W = W l by looking at 



P13P24 (be + dc)af 



P12P34 



(bf)(ae) 



1 + 



dc 

Ve 



1 + W; 
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PJ3P24 _ (EF)(AC + BD) = ! 

p' 12 ^ 34 (CF)(A£) +l J • 



The condition Eli = £12 l eac j s to the following: 

P24 P24 

- = -. ^ = tt/^h => wl = wi(i + wo. 

The rest of the relations (|35]) . see Figure [2T] are obtained by looking at 22a, 2a* 5 anc l 214. 

Conversely, transformation formulas (|38p imply that the vectors (pij) and (pL) in CP 5 coincide. 

Furthermore, we claim that (|4ip implies (|40p . The subtlety is that the modified partition 

function involves the function a. Let us set, see the right picture on Figure | 



Qi = P3 ~ Pi,Q2 =P2~ P3, <73 = P4 ~ P2, Qi = Pi ~ Pi- 
Then the denominators in (|40p for the inner graphs are 

M™ = a ^-^¥ >2 - ps c p '^ Pl (F i - p2 e pl - pi f P2 - pi = (cf) qi (bf) q2 (ad) q3 (ae) q4 , 

M /in = A Pl ~ Pi B P4 ~ P2 C P2 ~ P3 D P3 ~ Pl E P3 ~ P4 F P4 ~ P:i = (DE) qi (CF) q2 (BF) q3 (AE) q4 
Therefore, using gTJ, M™ = C qi+q2+q3+q4 M'™ = CM'™ since q x + q 2 + q 3 + <? 4 = 1. So, using 



411) again, we get (jlOjl . D 



o 

Figure 23: Shrinking a 2-valent vertex. 

There is another transformation, see Figure [23"1 which obviously does not change the face 
weights and the Poisson structure. Its inverse involves splitting of a high valency vertices by 
adding valency two vertices. 

4.3 Quantum integrable system 

Recall the quantum torus algebra Ta assigned to a lattice A with a skew-symmetric integral 
form (*, *), see Section [4.1.31 Each vector v € A gives rise to an element X v € Ta- 

Each Hamiltonian H a ^ a is a sum of different monomials with coefficients +1. So it can 
be upgraded uniquely to an element of the quantum torus algebra Ta assigned to a minimal 
bipartite surface graph T: 

H a ,:=^I„eT A (44) 

V 

where the sum is over the vectors v parametrising the monomials in the classical Hamiltonian 
H a)a . The element (|44p is the quantum Hamiltonian. Similarly one defines quantum Casimirs 
corresponding to the classical Casimir monomials. 
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Proposition 4.9. Let Y — >■ V be a spider move. Then the corresponding quantum cluster 
transformation preserves the quantum Hamiltonians and Casimirs. 

Proof. The Casimirs lie in the center and thus commute with ty q (X v ). 

For the quantum Hamiltonians we use a version of the strategy employed at Section 14. 21 
First, let us do a calculation at q = 1. Denote by qijPij (respectively qLpL) the part of the 
modified partition function involving the dimers connecting the vertices i and j to the outside 
vertices. Here pij (respectively j/ •) is the factor given by the matchings supported on the edges 
a, b, c, d, e, f (respectively, A, B, C, D, E, F) on Figure [22j 

We write A ~ B when A/B is independent of 

W,W 1 ,W 2 ,W 3 ,W i , and W, W{, W£, W£, W{. (45) 

Since the product of all face weights is equal to one, we have 

WW 1 W 2 W z W i ~ 1, W'W[W' 2 W Z W[ ~ 1. 

Below we work modulo these relations. The 1-cycle assigned to any non trivial Laurent polyno- 
mial in W, Wi, W 2 , W 3 , Wi has a non-trivial contribution of the edges a, b, c, d,e, f. 

The qij are Laurent polynomials in the face weights for T, considered modulo a common 
Casimir used to correct the partition function. They do not depend on the W's listed in (|45p . 
So we can calculate the contribution of these W's to the ratios of p's by looking at (I42p . For 

example, we get 

Pis be + dc e dc 

— ~ -^7- = 7 + T7 ~ w 3 + ww 3 . 

P12 bf f bf 

Indeed, the contribution of a, b, c, d, e, / on T to the 1-cycle assigned to W 3 is [e] — [/], and the 
contribution to WW 3 is [d] + [c] — [b] — [/]. These Laurent polynomials in W are uniquely defined 
modulo WW1W2W3W4 ~ 1. 

The same argumentation in general delivers the following result: 

— ~ W 3 + WW 3 , — ~ WW3W4, — ~ WW 2 W 3 , 
P12 P12 P12 



Similarly, we get 



P12 P12 

Pl2 P12 Pl2 

P12 P12 

Now it is straitforward to chack that the mutation at W formula amounts to mapping the five 
Laurent polynomilas in W's to the similar Lauren polynomials in W's. Moreover, the same is 
true for the quantum mutation maps, assuming that the Laurent polynomials above are written 
with the correcting ^-factors provided by the formula 



W3 1— ► W 3 (l + qW). q~ X W' z W[ t— ► W3WW4, . . . 
It remains to notice that WWiW 2 W 3 W4 1— > W'W[W' 2 W' z Wl U 
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Theorem 4.10. Let T be a minimal bipartite graph on a torus T. Then the quantum Hamilto- 
nians |^| ) commute. 

The quantum Casimirs generate the center. The quantum Hamiltonians H Q;a together with 
the quantum Casimirs provide a quantum cluster integrable system in O q {Xp{). 

Proof. One has 

[X v ,X w ] = (q^-qM)X v+w . 

Using the notation of Lemma 13.81 we have 

[X^X^] + [X^Xfr] = (q<^>l») - q (^))X, 1+ll2 + (q&w) - q^~^)X~^ 2 . 

By construction we have [i\ + ^2 = Mi + JU2- So the commutator is zero by Lemma 13.81 There- 
fore the quantum Hamiltonians commute. The claim about Casimirs is obvious since they are 
monomials. The last claim follows Proposition 14.91 and Theorem 13. 7l i). 

□ 

5 Resistor network model and its cluster nature 

5.1 Resistor network model 

Let G be a (not necessarily bipartite) surface graph on a torus T. Let c be a function on the 
edges of G with values in C*, called the conductance function, considered up to a multiplication 
by a non-zero scalar HI A graph with a conductance function on its edges is called a resistor 
network. 

5.1.1 Laplacians 

A line bundle with connection V on G gives rise to a vector space 

given by the sum of the fibers V v of the line bundle over the vertices v of G, and a Laplace 
operator A : V — > V defined by 

HfKv):=J^c(v,v')(f(v)-i v , v f(v')) 

where the sum is over neighbors v' of v, and i v 'vfW) ls ^ ne parallel transport of the vector f(v') 
to the fiber over v. 

Suppose now that the bundle with connection V on G is flat, that is, the monodromies 
around the faces of G are trivial. The connection still has non-trivial monodromies Z\, z<i around 
generators of iJi(T,Z). The determinant of A is a Laurent polynomial det A = P{z\,z?). It is 
symmetric, since A is Hermitian when \z\\ = \z2\ = 1 and c is a positive real function: 

P(z 1 ,z 2 ) = P(z^\z 2 1 ). (46) 

The Laplacian determinant is a weighted sum of combinatorial objects called cycle-rooted 
spanning trees (CRSFs), see JK]: these are maximal subsets of edges of G in which each connected 



7 Usually in applications the conductance function is a positive real- valued function. 
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component contains a unique cycle, and this cycle has nontrivial homology on the torus. The 
weight of a CRSF is the product over its edges of the conductances, times the product over its 
cycles of 2 — w — 1/w where w is the monodromy of the flat connection on that cycle. 

Theorem 5.1 ([K]). We have 

detA= Yl II c ( e ) II 2-w(rf)-l/w(i). 

CRSFs C eeC cycles 7 

Two torus graphs G\ , G2 with conductances c\ , C2 are said to be electrically equivalent if 
their associated Laplacian determinants P\{z\,z-i) and Pi{z\ y z-i) are proportional (that is, their 
ratio is a monomial in z\,z% ). A weaker notion is topological equivalence: torus graphs G\,Gi 
are topologically equivalent if for (any choice of) positive conductances, the Newton polygons of 
Pi(zi,Z2) and Pi{z\i&l) are the same. 

5.1.2 Dual graph 

Given a resistor network on a surface, the dual graph G' is also a resistor network, where we 
assign conductance 1/c to an edge dual to an edge of conductance c. The Laplacian determinants 
on G and G' are related: 

Lemma 5.2. For a graph G with dual G' having the same flat line bundle, we have 

] [ c(e) J det A G / = det A G . (47) 

edges J 

In particular a resistor network on a torus is electrically equivalent to its dual network. 

Proof. Use Theorem 15.11 and note that, given a CRSF C on G, there is a dual CRSF C' on G' 
consisting of the set of edges of G' which are not crossed by an edge of C. It has the same 
homology class as C, and thus the same weight Ilcycics 2 — w — l/u>, and the contribution from 
the edge weights is IleeE c ( e )~ 1 times the product of the edge weights of C. □ 

5.1.3 Zig-zag paths for a graph G on a torus. 

A zig-zag path on a surface graph G is a path which turns maximally left or right at each vertex, 
and the left/right turns alternate. A zig-zag path on G is unoriented. The collection {a} of 
zig-zag paths assigned to a graph G on a torus has the following property: 

Lemma 5.3. Each zig-zag path intersects an even number of zig-zag paths (counted with mul- 
tiplicity). 

Proof. Each vertex of a zig-zag path is either on the left or right of the corresponding path in 
the medial graph, so a zig-zag path has an even number of steps. At each step it crosses a single 
other zig-zag path. □ 

Conversely, any collection of smooth unoriented loops {a} on a torus, in general position, 
satisfying Lemma [5731 and considered up to isotopy, arises as the set of medial strands of a graph 
on the torus. Namely, one can color the complementary domains black and white so that no 
domains of the same color share an edge. We declare the black domains to be the vertices of 
G, and connect two vertices of G by an edge if the corresponding domains are separated by an 
intersection point of two zig-zag paths. The white domains are the faces of G. 
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5.1.4 Elementary transformations. 

Elementary transformations are local rearrangements of a resistor network which preserve elec- 
trical equivalence. There are four types: Y — A transformations, parallel and series reductions, 
and removal of dead branches. 

Recall that a Y- /^-transformation of a resistor network is given by replacing a triangle A 
with conductances a,/3,7 by the graph Y with conductances a,b,c, see Figure l24"l related as 




Figure 24: An example of a Y-A-transformation. 



follows: 



be 



a 



a 



a + b + c 

aj3 + /?7 + 7a 



, /3 



or 



06 



a + b + c 



, a/3 + £7 + 7a 

= , c 



a + 6 + c 

a/3 + /?7 + 7a 



(48) 
(49) 



a /3 7 

The graph transformation induced by moving a zig-zag strand across a crossing of two other 
zig-zag strands is a Y-A move. 

Two more elementary transformations of graphs with conductances, called parallel and series 
reductions, are shown in Figure! 



C 2 " 




c 1+ c 2 


c l 

> 
c 2 




c l c 2 




c 1+ c 2 





Figure 25: Parallel and series reductions. 

A last elementary transformation is removing a dead branch, which consist of contracting to 
a vertex any subgraph which is attached to the rest of the graph by a single vertex. 

A torus graph is said to be minimal if the lifts of any zig-zag path to the plane (the universal 
cover of the torus) does not intersect itself, and the lift of two zig-zag paths to the plane intersect 
at most once. 

Theorem 5.4. Every torus graph G with conductances c can be reduced to a minimal torus 
graph by a sequence of elementary transformations. 

Two minimal torus graphs G\,G2 are topologically equivalent if and only if G<i is related to 
G\ or to its dual graph G'i by Y-A-moves. 



We prove this at the end of Section 15.2.41 The corresponding statement for electrical equiv- 
alence is false: there are minimal graphs which are electrically equivalent but not related by 
Y — A moves. In fact the electrical equivalence class is a subvariety of the Hamiltonian dynam- 
ical system; see below. 
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5.2 The associated dimer model on a torus 

5.2.1 The bipartite graph assigned to a graph on a torus. 

Drawing the graph G and its dual graph G' simultaneously one can make a bipartite graph Tq 
on the torus whose white vertices are the edges of G, black vertices are the vertices and the faces 
of G, with an edge for each adjacency, see Figure [2j The graph Tq has the same number of white 
and black vertices since Euler's formula on a torus reads E = V + F. One can show |KPW| that 
the graph Tq admits dimer covers: in fact these are in bijection with CRSFs in which both the 
cycles and dual cycles (the dual of a CRSF is also a CRSF) have a chosen orientation. 

5.2.2 The resistor network subvariety. 

Let c be a conductance function on the edges of G, considered up to a non-zero scalar. It gives 
rise to a line bundle with connection V(c) on Tq, defined as follows. The fibers of V(c) at the 
vertices of Tq are identified with C. So the connection is defined via the edge weights assigned 
to the oriented edges of Tq- The weight of an edge incident to a vertex corresponding to a face 
of G is 1. The weight of an edge going from a vertex of G to a midpoint of an edge E of G is 
the conductance ce assigned to that edge. 

The monodromies around a pair of generators of Hi(Y, Z) provide an isomorphism of the 
moduli space of line bundles with flat connections on Tq with (C*) 2 . Taking the tensor product 
of V(c) with a line bundle with a flat connection on Tq we get a line bundle with connection 
on Tq- We denote by IZq the moduli space of line bundles with connections on the graph Tq 
constructed in this way, and call it the resistor network subvariety. It is a closed subvariety 

1Z G C £r G , dim7e G = E G + 1 

where Eq is the number of edges of G (note that multiplying all conductances on G by a constant 
gives a gauge equivalent connection on Vq). 

5.2.3 Determinant of the Laplacian and the dimer model partition function 

Using the conductances on G to define the connection on Tq, we have the following relation 
between the dimer partition function and det A. 

Lemma 5.5. The partition function for the dimer model on Tq is equal to detAg (up to a 
monomial factor in Z\,z-i)- 

Proof. Let G' be the dual graph to G with conductance 1/c on an edge dual to an edge of 
conductance c. Recall Lemma [4"T1 relating the Laplacian determinants on G and G' . 

These two operators A^, Ac acting on the black vertices of Tq are related to the Kasteleyn 
matrix K as follows: 

where K\,K2 are two copies of the black-to-white Kasteleyn matrix with appropriate gauge, 
defined as follows. 

Embed G in a Euclidean torus with straight edges; this allows us to associate to each directed 
edge 6162 a unit-modulus complex number C(^i 3 ^2) representing its direction. Letting w be the 
vertex of Tq in the center of this edge we then define K\{bi,w) = c(bi,b2)C(°i,°2) where c 
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is the conductance of edge 6162. If 63 is a black vertex of G' in a face of G adjacent to edge 
6162 then Ki(b3,w) is 1 times the unit modulus complex number in the direction from 63 to w 
and perpendicular to ((bi, 62). One can now check that the monodromy around a face of T^ is 
negative real (for positive conductances) and thus K\ is a Kasteleyn matrix for Tq. One defines 
K% similarly but for the dual graph G'. Now (|50p follows by inspection of the various cases: 
—K^Ki^bjb) is the sum of the conductances of edges emanating from b; —K^K^b^b') = if 
b G G and b' £ G' or vice versa; —K^Kiibx, 62) = — c if 61, 62 are adjacent vertices of G and the 
edge connecting them has conductance c; similarly for vertices adjacent in G'. 

Finally, det K2 , det K\ are both the dimer partition function (up to constants which depend 
on our choice of gauge) . Thus from (|47p and (|50p we have det A = det K\ up to multiplicative 
factor which is a monomial in z\, z%. □ 

It follows from ()46|) that the Newton polygon of the dimer model on Tq is centrally symmetric 
with respect to the origin. 

5.2.4 Zig-zag paths and the Newton polygon 

Zig-zag paths are defined for both a graph G on a torus and the associated bipartite graph Tq', 
zig-zag paths on G are unoriented. Zig-zag paths for Tq are oriented. A zig-zag path a on 
G corresponds to a pair of zig-zag paths on Tq, isotopic to a and equipped with the opposite 
orientations. In this way we get all zig-zag paths on Tc- 

The geometric Newton polygon N(G) is described as follows. Each zig-zag path a on G gives 
rise to a pair of homology classes ±[a] € H\(T, Z), where [a] is the homology class of the path a 
equipped with an orientation. There is a unique, centrally symmetric with respect to the origin, 
convex integral polygon 

N(G) eiTi(T,Z) 

with the sides given by the vectors ±[a] for all zig-zag paths a on G. Indeed, there is evidently 
a unique up to a translation convex integral polygon in Hi(T, Z) with these sides, which is 
centrally symmetric. Thanks to Lemma 15.31 it can be centered at the origin. 

Clearly the polygon N(G) coincides with the geometric Newton polygon N(Tg) assigned in 
a similar way, using the homology classes of zig-zag paths, to the bipartite graph Tq- Thus 
Lemma 15.51 implies 



Lemma 5.6. Let G be a minimal graph on a torus. Then the geometric Newton polygon N(G) 
coincides with the Newton polygon of det Aq- 

Proof of Theorem 15.41 By Lemma 15.61 t wo equivalent minimal graphs G and G' have the 
same geometric Newton polygons, whose boundary sides determine the homology classes of the 
zig-zag strands. Thus the collections of their zig-zag strands are isotopic. The isotopy is realized 
by moving zig-zag strands through zig-zag strand crossings, and thus the corresponding graphs 
are related by a sequence of Y-A moves. 

To reduce a graph G to a minimal graph, consider the zig-zag strands for G, and isotope 
each to a straight geodesic through a generic isotopy. Each type of singularity which resolves 
corresponds to an elementary transformation: if it is a parallel, series, or dead branch move, it 
may change the topology of the set of strands, but also strictly reduces the number of crossings 
of strands. Eventually one arrives at a graph with a minimal number of strand crossings and 
the remaining moves are Y — A moves. □ 

44 



5.2.5 Equations denning the resistor network subvariety. 

Denote by Vq, Eq, Fq the number of vertices, edges and faces of the graph G, and similarly 
Vr, Ep, Fr those for Tq. Then Fp = IEq. Indeed, every edge E of G is incident to four faces 
of Tq, and each face of Tg is counted twice. So one has 

dim£ r = F r + 1 = 2E G + 1. 

So the codimension of the resistor network subvariety IZg is Eq. On the other hand, a line bundle 
with connection on Tq which belongs to the subvariety IZg satisfies the following conditions: 

1. The monodromies around the faces of the graph G and the graph G' are equal to 1. 

2. Given a loop a on the graph G and a loop a' on the graph G' whose homology classes in 
.Hi(T,Z) coincide, the monodromies around these loops also coincide. 

Remarks. The first condition tells that for each black vertex b on Tq, the product of the face 
weights for the faces adjacent to b is equal to 1. 

Secondly, the second condition implies the first: to see this, alter a loop a by adding to it a 
loop going around a face of the graph G, and similarly for G' . 

Finally given the first condition, the second condition imposes only two more equations on 
the connection: one for each generator of H\{Y, Z), since the first condition already implies that 
the monodromies around homologous paths, both in G or both in G', are equal. 

Since the product of the monodromies over all faces of G, as well as over all faces of G' , is 
1, we get Eq equations for the subvariety TZq C Ct g : 

(V G + F G -2) + 2 = E G . 

These conditions characterize the subvariety TZq in £p ; 

Lemma 5.7. A line bundle with connection on Tq arises from a resistor network on G if and 
only if it satisfies the conditions 1 and 2 above. 

Proof. We have seen already that the line bundle with connection on Tq arising from a resistor 
network satisfies the conditions 1 and 2. 

Let us prove the converse. Notice that the condition 1 says that we have a flat line bundle 
on the torus T assigned to the graph G, and another flat line bundle on the torus assigned to 
the graph G'. The condition 2 then says that these two flat line bundles are isomorphic. So, 
untwisting by a flat line bundle on the torus, we can assume that our line bundle has trivial 
monodromy around any loop on G and around any loop on G'. 

Let us trivialize the line bundle on Tq on the dual graph G'. Take an edge Eq of G and 
assign it conductance 1. Recall that the medial graph of G is the graph with a vertex for every 
edge of G and an edge whenever two edges of G share a vertex of G. The parallel transports 
on Tq determine a curl-free flow on the medial graph of G; the value of the conductance of an 
edge E is then the integral of this flow along a path in the medial graph from Eq to E. So we 
have reconstructed conductances up to a non-zero common factor. □ 
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5.2.6 The resistor network subvariety is Lagrangian. 

An irreducible subvariety of a Poisson variety is Lagrangian if the generic part of its intersection 
with the generic symplectic leaves are Lagrangian subvarieties of the symplectic leaves. 

Theorem 5.8. IZq is a Lagrangian subvariety of the Poisson variety Cv G - 

Proof. The following crucial lemma implies that the intersection of IZq with a generic symplectic 
leaf is coisotropic (that is, the symplectic form is zero on this intersection). 

Lemma 5.9. The Poisson brackets between the defining equations 1-2 of the resistor network 
subvariety IZq are equal to zero. 

Proof. Let us show first that the monodromies around the black vertices of the graph Tq Poisson 
commute. Let 61,62 be two black vertices of Tq corresponding to the faces F\,F% of G. Let 
Wb 1 ,Wb 2 be the corresponding monodromies. If the faces F\,F2 are not adjacent, the Poisson 
bracket is obviously zero. If they are adjacent, they share an edge E, which has, by construc- 
tion, two black vertices at the ends and one white vertex at the middle, see Figure El Thus 
{Wi^, W& 2 } = 0. The argument in the dual case when 61,62 are two black vertices of Tq corre- 
sponding to the vertices of G is the same. Finally, let 6 be a black vertex assigned to a face F 
of G, and b' a black vertex assigned to a vertex v of G. Then the face F of G and the face F' 
of the dual graph G' corresponding to v either do not intersect, or intersect at two points. See 
Figure [2] where the boundary of F is a black loop and the one for F' is a dashed loop. In the 
former case the Poisson bracket is evidently zero. In the latter case, the contributions of the 
intersection points come with the opposite signs, and thus cancel. □ 

Continuation of proof of Theorem \5.8[ Let 2m be the number of Casimirs for the Poisson variety 
£r G - It equals to the number of primitive boundary intervals of the Newton polygon N(Tg) 
and thus is an even number, since the polygon is centrally symmetric. Since the product of all 
Casimirs is 1, the dimension of the generic symplectic leaf on Cy g is 

2E G + 1 - (2m - 1) = 2{E G - m + 1). 

The Casimirs are monodromies along the zig-zag paths, and the latter on the graph Tq come in 
pairs: {zi,z[}, ..., {z m ,z' m }, so that {zk,z' k } have the opposite homology classes, and thus the 
corresponding Casimirs Ck,C' k are inverse to each other, providing a relation CkC' k = 1. Since 
the product of all Casimirs is 1, the last relation C m C' m = 1 follows from the others, while the 
first m — 1 are independent. 

Since dimlZ equals to m — 1 plus half of the dimension of the generic symplectic leaf, it is 
sufficient to show that the generic symplectic leaves intersected by 1Z are parametrized by (no 
less than) m — 1 parameters. Let us assume that there is a relation C" 1 . . . C^ m = 1. Let us 
show that m = . . . = n m . Take a pair of Casimirs C a and Cb corresponding to non-parallel sides 
of the Newton polygon. Then the corresponding zig-zag strands intersect. So there is an edge 
E shared by them. Altering the conductance at the edge E we alter just the Casimirs C a and 
Cb- This implies that n a = n^. Therefore m = . . . = n m . □ 
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5.2.7 Quantized cluster monomial Lagrangian subvarieties 

The subvariety 1Z is a cluster monomial Lagrangian subvariety - it is a Lagrangian subvariety of 
a cluster Poisson variety X , which is given in any cluster coordinate systems by equating certain 
monomials to 1. 

According to the correspondence principle of quantum mechanics, a symplectic manifold X 
should give rise to a Hilbert space Hx - the quantization of X. Furthermore, a Lagrangian 
subvariety L C X should provide a line in Hx- Here is how the story goes in the cluster 
situation. 

The quantization of a cluster Poisson variety X |FG2j provides a collection of Gelfand triples, 
parametrized by the eigenvalues \ °f the Casimirs: 

Here H x is the Hilbert space discussed in Section 11.2.21 The modular double algebra, given by 
the tensor product of quantum algebras at h and 1/h: 

O q {X) ® O q y (X), q = e™\ q v = e™'\ (51) 

is realized by unbounded selfadjoint operators in the Hilbert space H x . The subspace S x is the 
maximal subspace of H x where the algebra ([5T]) acts, equipped with a natural topology. The 
space S* is its topological dual. One can show using the methods of loc. cit. that any cluster 

A. 

monomial Lagrangian subvariety L <Z X gives rise to a functional 

defined uniquely up to a scalar by the condition that tpi, is annihilated by the quantized monomial 
equations in the algebra O q {X), and their counterparts in the algebra O q v(X). 

Here is the basic example. Let H = L 2 (M) be the Hilbert space of functions f(x). Consider 
unbounded operators 

Xf(x) := f(x + inH), Yf(x) := e x f(x), 

X^f(x) := f( x + in), y v /(x) := e x / h f{x). 

Then 

XY = q 2 YX, X V Y V = q v2 Y v X v . 

So we get a representation of the modular double of the quantum torus algebra with two gen- 
erators. The subspace S C % consists of functions f(x) with exponential decay at x — > ±oo, 
whose Fourier transform has the same property. By the Pay ley- Wiener theorem f(x) has an 
analytic continuation to C 

The equation X = 1 defines a Lagrangian subvariety L in the torus (C*) 2 with coordinates 
(X, Y). The quantized equations are (X — l)f(x + iy) = 0, (X v — l)f(x + iy) = 0. For irrational 
h, they imply that f(x + iy) = f(x). This implies that the functional ip is a constant. 

5.3 Cluster nature of the resistor network space 

Recall that all our spaces are equipped with a free action of the two-dimensional torus T, 
and relating them to cluster varieties we have to consider their quotients by this action, or, 
equivalently, consider the action of all transformations on the face weights only. 

In Section 15.3.11 we use ' in the notation X' G and 1Z' G to indicate that, just like in Section 
14.1.61 we work with the quotients Xq/T and 7^-g/T. Moreover we do not define the spaces Ag 
and Bq, only the "quotients" A' G and B' G . 
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5.3.1 The basic commutative diagram for a graph G 

Given a graph G on a torus, let us define four split algebraic tori (that is, varieties (C*) fc for 
some k) A' G , B' G , 1Z' G , X g , fitting into a commutative diagram 



(52) 



The tori A' G , B' G , TZ' G , X' G are defined as the tori with the following canonical coordinates: 

• A' G : coordinates {^4/} at the faces / of the graph Tq. 

• B' G : coordinates {B b } at the black vertices b of the graph Tq. 

• 1Z' G : coordinates {C w } at the white vertices w of Tq, defined up to a scalar. 

• X' G : coordinates {Xf} at the faces / of Tq, whose product is 1. 

The C^-coordinates are nothing else but the conductances assigned to the edges of the graph 
G, which are identified with the white vertices of Tq. 

We define the maps between these tori by their action on the coordinates. Below ~ means 
an incidence relation. So b ~ / means that a black vertex b is incident to a face /, etc. 

First, the map p is the standard map [FGl^ 



p*X f :=l[A £ /°, (53) 



where the product is over faces g of Tq, and £f g is the bilinear form on the faces defined in 
Section 2. 
Next, 

j*A f :=Y[B b = B b B b ,. (54) 

6~/ 

Here b, b' are the two black vertices of the face / of the graph Tq. 
Then, we set 

«**/ ■= ^T- (55) 

Here (61,621^3)^4) are the four black vertices which lie outside of the face and incident to the 
two white vertices of the face, as shown on Figure 1261 

Next, 

i*X f := £*, (56) 

Here w, w' are the two white vertices of the face /, and w, b, w' go counterclockwise around the 
face, and the black vertex b is a vertex of /, Figure l27l 
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Figure 26: The map s*Xf. 
b, 




Figure 27: The map i*Xf := C w /C w >. 



Finally, 



„*r< Cl CJ 

q w '~ R R 

D C2 D CA 



(57) 



where c\ , c% , cs , C4 are the four black vertices incident to the white vertex w as shown on the 
Figure [281 

Define [FGj a 2-form on the space A' G by 

n A := ^2e fg dlogA f AdlogA g . 
f,g 

Proposition 5.10. i) The diagram \5Sfy is commutative. 

ii) The map j is an embedding. Its image is isotropic for the 2-form fi^. 
Hi) The subvariety i(jV G ) is defined by monomial equations. It is Lagrangian. 
iv) The subvariety s(B' G ) is Lagrangian in the symplectic variety 

U' G :=p(A G )cX G . 

Proof, i) One checks this by looking at the Figure [2HJ I27T [2"B1 
ii) Follows easily from the definitions, 
iii) The subvariety i(TZ G ) C X' G is given by monomial equations 

n x f =i - ( 5g ) 

/: Yldf=a-I3 
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Figure 28: The map q*C u 



Here equations are assigned to pairs of loops a, (3 on G which are homologous on the torus. The 
product is over a set of faces / of T G whose total boundary is a — j3. In particular, taking the 
empty loop /3, we get equations 



I I Xf = 1 for each black vertex b of the graph Fq- 
/~6 



(59) 



Indeed, it is easy to see that i* (|58j) = 1. Recall that E G is the number of the edges of G. It is 
also the number of the black as well as the white vertices of Tq- One has 

dim i(K' G ) = E G - 1. 

Indeed, there are E G independent equations: E G — 2 are given by (I59p . and two more by /3's 
generating Hi(T,7*). Since dim 1Z' G = E G — 1, we proved the first claim. The second is similar 
to the proof of Theorem 15.81 

iv) It is similar to the proof of Theorem 15.81 

□ 




B 6 m c 3 B < 




B 3 B 6 9 B 7 > Gi • B 



Figure 29: Transformation law for B's under a Y-A-move. 



5.3.2 Action of the Y-A-move on coordinates. 

A Y-A-move acts on the conductances coordinates Cj as in (148p - (|49p . Referring to Figure l29l 
we have 



B1B7 



B2B3 B 3 B! , B a B! BjB 2 , B 1 B 2 B 2 B S 

' R.. It, I!,. «.. ~T 



h£>7 _ ~ _ C\C 2 + C 2 C 3 + C3C1 b B 4 B B 5 + B B 5 B B 6 + B B 6 B a B 4 
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from which one obtains 



B 7 



B\B± + B 2 B 5 + B 3 B 6 



(60) 



The same equation is obtained starting from the other two conductances Ci,C-$. This equation 
is the so-called cube recurrence relation, see Section 16.31 below. 



Lemma 5.11. The Y — A move G — > G' is a composition of four cluster mutations. Precisely, 
there are commutative diagrams 



A' 



G 



B> 



G 



Xq 



n' 



G 



IfJ-A 



I HX 



G' 3 G' 



£V/-1I \ l\>f~il 



where the down arrows are birational maps induced by the Y — A move action on the B- and 
C-variables, given by i60\) and |^<g|)-|7P|), and the right ones are compositions of four cluster A- 
and X-mutations. 

Proof, i) See Figure [30j After this sequence of 4 cluster transformations, a little algebra using 




Figure 30: The Y — A move from spider moves. 
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Theorem 14.71 shows that the new X variables are: 



X 
X. 

X- 
X 
x \ 

x { 

X- 

X, 
X, 



(-1) 



(4) 



(4) 



(4) 
i 

(4) 



(4) 
6 

(4) 

7 

(4) 
3 

(4) 



1 + X x + XjX 3 
X 3 (1 + X 2 +X 1 X 2 ) 

1 + x 2 + x x x 2 

Xi(l + X 3 + X 2 X 3 ) 
1 + X 3 + X 2 X 3 

^(l+Xi+XiXg) 

X 4 (l + Xi+XiX 3 ) 
X1X3 
5 i + x 1 + x 1 x 3 

X 6 (1 + X 3 + X 2 X 3 ) 
X 2 X 3 

A7 

1 + X3 + X 2 X 3 

X 8 (l + X 2 + XiX 2 ) 

X1X2 



X 



1 + x 2 + x ± x 2 



If XiX 2 X 3 = 1, we will also have X\ X2 X 3 =1. Moreover if the conductances on 
the edges of the triangle opposite faces 1,2,3 are Ci,C 2 ,C 3 respectively then X\ = C 2 /C 3 , 
X 2 = C 3 /Ci,X 3 = C\/C 2 . One checks that the new conductances are those satisfying (J4TJ|h 

ii) It follows by comparing formulas (|48p -(]49p. D 



It follows that the diagram (J52J) gives rise to a similar diagram of spaces, obtained by gluing 
the tori assigned to minimal equivalent graphs G via the corresponding birational maps given 
either by the Y-A-move or the corresponding sequence of four cluster mutations. 



6 Discrete cluster integrable systems 

6.1 Discrete cluster integrable systems on a square grid on a torus 

Let us consider the dimer model on a square grid on a torus. Such examples are constructed as 
follows. Take the standard square grid bipartite graph in Z 2 , and map it onto torus by taking 
the quotient by a rank two subgroup of the (bipartite-coloring preserving) translation group 
(that is, (1, 1)Z + (1, -1)Z), see Figure! 





% / J 



Figure 31: Inserting two white 2-valent vertices, and performing the spider move in the dotted 
region we get a version of the spider move on Figure [ 
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Let us construct an element A of the cluster modular group acting by a non-trivial automor- 
phism of the corresponding cluster Poisson variety Xn- 

There is a version of the spider move, see Figure 132"! given by a composition of two elementary 
moves on Figure [3T1 





Figure 32: A version of the spider move. 

The squares of the grid are of two types: one, which we refer to as the *-faces, have the 
black bottom left vertex. Let us do the Figure [321 version of the spider move at all *-faces. The 
result is illustrated at the right of Figure [Ml Then we shrink the 2-valent vertices, as shown on 
Figure 1231 We get again a square grid, but now the colors of the vertices flipped, see the bottom 
of Figure [Ml So performing the same sequence of moves — the Figure [321 version of the spider 
moves at the *-faces and shrinking the 2-valent vertices — we get the original grid. The resulting 
cluster transformation A is non-trivial. Thanks to Theorem 14.71 it commutes with the quantum 
Hamiltonians. 




Figure 33: Discrete cluster integrable system on a square grid on a torus. 

This transformation A is called urban renewal in [S]; its relation with the Hirota bilinear 
difference equation [Hir| is discussed there. See also [Kor| . 

The cluster transformation A together with the quantum Hamiltonians form a discrete quan- 
tum integrable system — we think of A as of an evolution with a discrete time. In the quantized 
picture [FG2] the operator A gives rise to a unitary operator in the Hilbert space of the system 
commuting with the quantum Hamiltonians. It is a composition of the quantum dilogarithm 
intertwiners parametrized by the faces. 
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Figure 34: Constructing an element A of the cluster modular group. 

6.2 Cluster automorphisms for resistor networks on a torus. 

One also finds cluster automorphisms for cluster varieties associated to resistor networks on a 
torus. The pair (X,H) assigned to an equivalence class of minimal graphs G on a torus has 
a large group A of cluster automorphisms, isomorphic to (a quotient of) Z s_2 , where 2s is the 
number of sides of the Newton polygon N{G). It is constructed as follows. (This material is 
also found in |HSJ.) 

For each pair E, E' of centrally symmetric sides of the Newton polygon N(G) of a minimal 
graph G on a torus there is a collection ai,...,atk of isotopic zig-zag strands on G which do not 
intersect due to minimality of G: their (well defined up to sign) homology class in -Hi(T, Z) is a 
primitive vector on E. We assume that they are written in their natural cyclic order provided 
by the embedding to the torus. Isotope them cyclically (ai,...,ak) — > («2j ••-,«&, «i), leaving 
the other zig-zag paths intact. Present this as a sequence of Y-A-moves corresponding to the 
zig-zag strands moving through simple crossings of other zig-zag strands, which transforms the 
graph G to itself. Take the product of the cluster transformations assigned by Lemma 15.111 to 
the Y-A-moves. We arrive at a cluster transformation, realized by a birational isomorphism of 
the pair {X,1Z). It is an element of the cluster modular group. The automorphisms assigned to 
different sides commute, and their product is 1. 



6.3 The cube recurrence discrete cluster integrable system 



Let us apply the construction of Section 16.21 to the honeycomb graph on a torus. In this case 
there are three families of parallel zig-zag loops. Moving cyclically one of the families we get 
a cluster automorphism. We show that its action on the S-variables is given by the so-called 
cube recurrence. 



54 



A function a : Z 3 — > C satisfies the cube recurrence if 

a x+l,y+l,z+l a x,y,z = a x+l,y,zO>x,y+l,z+l + 0,x,y+l,zdx+l,y,z+l + 0>x,y,z+l a 'X+l,y+l,z- (61) 

Note that if a is defined and nonzero for 0<x+y+z<2 then equation (|61[) defines it for 
x + y + z = 3, also for x + y + z = — 1, and by iteration for all (x, y, z) 6 Z 3 . 



(1,-1,1) (1,)W5 (M,0) 



(0,-1,2) (0,Ai^ (141) (1,0,0) (1,0-1) OA.-1) 



(-li-l^ (0,0,2) (0,0,1) (0,0,0) (1,1,0) (1,1.-1) (l,a-3) <2,2,-2) 



(-1,0,1) (0^1) (0,1,0) (0.L-1) (1,2,-1) 



(-1VM3) (0^0) 



Figure 35: Vertex and face labeling for the cube recurrence. 

Let G be the honeycomb graph, thought of as a resistor network. Vertices of G are the 
orthogonal projection to the plane {x + y + z = 0} of the set {(x, y, z) GZ 3 | x + y + z = or 2}; 
edges connect vertices whose coordinate sum is to those nearest neighbors whose coordinate 
sum is 2. We can likewise index faces of G by (x, y, z) satisfying x + y + z = 1. See Figure 1551 
Let Vk be the set of vertices or faces with x + y + z = k. 

Assign variables B X ^ V)Z to the vertices and faces of G. 

Perform a Y — A move at each vertex v x ^y tZ where x + y + z = 0. The graph becomes 
the regular triangulation with vertices at V2, faces at V\ and new faces at V3. Using (|60j) . the 
.B-value at a new face is obtained from (|6ip . See Figure [36l 

Next, perform A — Y moves on all faces in V\. The resulting graph is again a honeycomb, 
with vertices V2, V$ and faces at V3. Iterating this procedure we produce variables B XiV;Z for all 
(x, y, z) £ Z 3 and these variables satisfy the cube recurrence. 

Equivalently, we can look at zig-zag strands and move cyclically one of the families of parallel 
zig-zag strands, say the vertical ones on Figure [35l Then each strand will cross two different sets 
of crossings of other two familes, which can be thought of as the vertices of the grid, resulting 
into Y — A and A — Y moves described above. 

By Lemma 15.111 the composition B of these Y — A moves and A — Y moves is a cluster 
automorphism, i.e. an element of the cluster modular group, of the dimer model phase space 
assigned to the graph G. Therefore it gives rise to a discrete cluster integrable system, with 
the discrete time evolution given by B. Thanks to (|60p . the cube recurrence is obtained by 
restricting the cluster automorphism B to the isotropic subvariety B. 
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Figure 36: Steps of the cube recurrence. The black dots are the vertices v xy z with x + y + z = 0. 

7 Newton polygons, toric surfaces and integrable systems 



7.1 Newton polygons and toric surfaces 

Let us briefly recall this well known connection. See [F] for more details. Let T = (C*) 2 . Denote 
by X*(T) the group of its characters, isomorphic to Z 2 . 

A Newton polygon TV is a convex polygon in the real plane -X*(T)]r := X*(T) <8>M = IR 2 with 
vertices at the lattice X*(T). It gives rise to a polarized projective toric surface A/*, possibly 
singular. The polarization is given by a T-equivariant line bundle C on M. Shifting the Newton 
polygon by a lattice vector I £ X*(T) amounts to changing the equivariant structure of the 
line bundle C by multiplying it by the trivial line bundle with the torus T action given by the 
character xi assigned to I. 

Each edge E of TV gives rise to a ray pe in the dual space X*(T)r, consisting of the normals 
to the edge E looking inside of TV. These rays form a fan A^v- It subdivides the plane X*(T)r 
into a disjoint union of cones. Each cone a provides the dual cone <r v , and hence a semigroup 
R a := <t v fl X*(T) and an affine surface U a := Spec T\R a \. Gluing the surfaces U a so that 
U ai na 2 = U ffl fl U a2 we get the surface J\f. 

The torus T acts on J\f. There is a unique 2-dimensional orbit, Uq, where is the origin. 
Its complement is the divisor at infinity Moo, whose shape is described by the polygon N: the 
irreducible components De are projective lines parametrized by the edges E of N, and their 
intersections match the vertices of N. The kernel of the action of the torus T on the component 
De is the cocharacter Xe '■ ^* "" ^ T(C) assigned to the conormal line to the edge E. 

There is a canonical isomorphism, serving as a definition of the line bundle C: 



■IN 



C^O(D), D = Y j a E D E . 



E 



where the coefficient cle € Z is the distance from the origin to the edge -EO The torus T action 
on C comes from the natural T-action on O(D), provided by the action on J\f preserving the 
divisor D. Shifting the polygon A by a vector I £ X*(T) amounts to multiplication of the 
isomorphism ijy by the character xu i- e - ^n+i = iNXl- Using the isomorphism zjy, we have 



H°{Af,£) ^> H°(Af, O(D)) c C(T) = C{z l} z 2 ), 



(62) 



where {z\,z-i) is a basis in the group of characters of T. It is identified with the space of all 
Laurent polynomials with the Newton polygon N. 



The latter is the value at the origin of the unique surjective affine linear function oe '■ X t (T) 
vanishes on E and is non-negative on TV. 



Z which 
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So given a pair (A/", C) , the polygon N is recovered as the convex hull of the set of characters 
of the action of the torus T on H°(M, £). 

Denote by u the canonical up to a sign invariant 2-form on M — Moo- Given a basis (zi, z<i) 
of characters of the torus T, it can be written as 

u) = d log z\ A d log zi . 

So it is a 2-form on M with logarithmic singularities at Moo- 



7.2 The spectral data fKQ] 



A parametrization of a degree d divisor C = Yli=i n * c * by a set 5 is a map 

v : 5 — ► {ci, ..., Cfe} such that |zv _1 (q)| = n^. 
Let |£| be a linear system of curves on M given by the zeros of sections of C. 

Definition 7.1. A spectral data is a triple (C,S,v) where C is a genus g curve on the toric 
surface M from the linear system \C\, S is a degree g effective divisor on Co = C n No, and 
v = {^e} are parametrizations of the divisors De fl C. 

The moduli space Sn = S parametrizes the spectral data related to N . 

Let C be a curve on M from the linear system \C\. Its intersection with N — Moo is defined 
by a Laurent polynomial P(zx, z<i) with the Newton polygon N. The intersection of C n De is 
determined by the sum of monomials of P supported on the edge E. The ratio of monomials 
assigned to two integral points on the sides of Af are Casimirs. So the Casimirs are naturally 
assigned to the primitive boundary edges of N, and their product is equal to 1. 

Proposition 7.2. Given a black vertex v of a minimal bipartite graph T assigned to N , there 
is a natural rational map 

Kr, v ■ ^ — ► S. 

Proof. A point of X parametrizes a line bundle with connection on T, encoded by the Kasteleyn 
operator K. Its determinant is a section of the line bundle £ ® L. The divisor of zeros of the 
section is the spectral curve C on M . The black point b defines a section of a line bundle on C. 
Its degree is the genus g of C, and its divisor is the divisor S. Homology classes [ag] of zig-zag 
loops on r are in bijection with the sides of the Newton polygon N. There is a parametrization 

7s : {Zig-zag loops on T in the homology class [ag]} ^> De n C, (63) 

uniquely determined by the following condition. The function given by the monodromy around 
a zig-zag loop a coincides with the function given by the coordinate of the boundary point 

IE (a). D 

Denote by S x the subvariety of S given by the condition that C C\ Moo is a given divisor 
Coo. Using the parametrization (f63j) . the divisor Coo provides a subvariety X x in X defined 
by the condition that the monodromies around zig-zag paths are given by the coordinates of 
the corresponding boundary points. Since this condition uses only the boundary divisor of the 
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spectral curve, there is a similar subvariety B x in B. The map Kr> is evidently fibered over B x , 
providing a commutative diagram 






\ 



6, 



Theorem 7.3. T/ie image of the map ky,v is open. The map k^,v is a finite map over the 
generic pointy] 

See below for the proof. Since both spaces have the same dimension, the first claim implies 
the second. 

7.3 The real positive points [KQ] 

Since X is a union of cluster coordinate tori, and the transition functions between different 
cluster coordinate systems are subtraction free rational functions, there is a well defined set 
A'(R>o) of its real positive points. It parametrizes connections on a minimal bipartite graph T 
associated with N whose face weights are positive numbers, and does not depend on the choice 

of r. 

Harnack curves were defined classically in [Hj as real curves P(z±, z%) = whose real compo- 
nents satisfy certain topological constraints. A simple characterization of Harnack curves, due 
to |MR| . is that they are the curves which maximize the area of their amoeba (which is the 
image of the curve under the map C 2 — > M? given by (^1,^2) l— >■ (log | ^1 1 , log | ^2 1 ) ) among all 
curves with given Newton polygon. 

The results of Kenyon-Okounkov |KO] can be interpreted as follows: 

• The real positive part i3(M>o) := tt(X(M. > q)) parametrizes Harnack curves on J\f which 
lie in the linear system \C\. Such a genus g Harnack curve has g ovals Oi matching the 
integral points inside of the Newton polygon N . 

• The space Af(R>o) parametrizes the data {A Harnack curve C on J\f from \C\, a choice 
of a point on each oval Oi, a choice of parametrization of the intersection of C with each 
divisor at infinity}. 

Therefore the fibers of the fibration 

7T+ : *(M >0 ) — ► B(R >0 ) 

are real (/-dimensional tori. They are the Liouville tori of the integrable system restricted to 
t>o)- 



7.4 The Hamiltonians 

Corollary 7.4. Forgetting the parametrization, the space B projects to an open part of \C\. 

Proof. Thanks to Section ITTBI the image of B is an algebraic subvariety of \C\, and its real part 
contains an open subset of the real locus of \C\. □ 



In the sequel to this paper we prove that this map is a birational isomorphism. 
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Corollary 7.5. The Hamiltonians are independent on X . 

Proof. Denote by \C\ X the moduli space of curves from the linear system \C\ with a given 
intersection C^ with Afc, . Thanks to Corollary 17.41 B x is a finite cover over on open part of 
\C\ X . We claim that the Hamiltonians are local coordinate functions on \C\ X . Therefore they 
are independent on X. 

Let us prove the claim. Identify the sections of C with Laurent polynomials in (zi, z-i) with 
the Newton polygon N. Let V^ be the subspace of Laurent polynomials in (z\, z-i) whose Newton 
polygon is the interior of N. Adding such a polynomial to a section of C, we do not change the 
intersection with A/"oo of its divisor of zeros. So we get an action of the vector space V% on \C\ X . 
This action is transitive. The Hamiltonians are given by the coefficients of the monomials in 
V%. So their differentials at a point C G \j0\ x form a basis of Tq\C\ x . □ 

Proof of Theorem 17.31 The restriction of the map Kr> to X x (M. > q) is an isomorphism. Thus 
its image is open in S x . □ 

8 Appendix: Calculating the Poisson structure on £p 

Let us give a local description of the intersection pairing (J3j> as a sum of local contributions 
assigned to the vertices of T. 

8.1 The local pairing. 

The star of a vertex v in a ribbon graph, Figure [37] is a vertex v together with the cyclically 
ordered edges incident to v. 




Figure 37: The star of a vertex on an oriented surface. 

The star of a vertex v gives rise to an abelian group A„, given by Z- linear combinations of 
the oriented edges Ei at v of total degree zero: 

A v = {^2 n i Ei I y^rij = 0, m € Z}, 

— >• 
where Ei is the edge Ei oriented out of v. Clearly A^ is a torsion free abelian group of rank 

-> — > 
m — 1, where m is the valency of v. It is generated by the "oriented paths" Ei — Ej, where 

— > 
— Ej is the edge Ej oriented towards v. 

Lemma 8.1. There exists a unique skew symmetric bilinear form 

5 V : A v A A„ — > -Z 
such that for any triple Eq,E\, E<i of edges as in Figure [321 on ^ has 

5v{lr A 7z) = £> Tr -=Ei - E , li :=E 2 - E . (64) 
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*1 



Y 

Figure 38: One has 5 v (j r ,ji) = \. 

Proof. Let {xi, ...,x m } be the cyclically ordered set of the endpoints of the edges sharing the 
vertex v. Let 7$ be the oriented path XiVXi+i. Then the elements 71, ...,7m generate A„ and 
satisfy the only relation £V ji = 0. We define the local pairing by setting 

SvijiAji+i) = -. 

Then ^2 i ji is in the kernel, so the local pairing is well defined. It is easy to see that it satisfies 
formula flM}. □ 

8.2 The global pairing 

Let L\, Li be two oriented paths on T. Let v be a vertex shared by them. Each path Li gives 
rise to an oriented path at the star of v, providing an element l{ £ A v . Let e v (L\,L2) be the 
value of the local pairing 8 V on l\ A l<i. 

Definition 8.2. Let L\,L2 be two oriented loops on a bipartite surface graph T. Then 

e(L 1 ,L 2 ) := y^sgn(v)e v (Li, L 2 ), 

V 

where the sum is over all vertices v shared by L\ and L2, and sgn(u) = 1 for the white vertex v, 
and —1 for the black vertex. 

^2 \ /Y 4 



^1 ^^3 

Figure 39: One has 5(71,72) = —1, 5(73,74) = 0. 

Proposition 8.3. The bilinear form e(L\,L2) from Definition {KM coincides with the bilinear 
form defined using the conjugated graph T in Section \l.l.l\ 

Proof. It is sufficient to show that in the examples on Figures [39] and 00] the sum of the local 
contributions coincides with the local intersection number of the corresponding paths on the 
conjugated graph. The latter clearly gives the results presented on Figures [39] and [40] So one 
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Figure 40: One has £(71,72) = 1, £(73,74) = 0. 

needs to show that the former gives the same result. For example, to check that £(71,72) 
for the left picture on Figure 2UJ notice that on Figure 01] we have 
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Figure 41: Proving £(71,72) = 1. 

Example. The face orientation provides the counterclockwise orientation of the boundary 
of the face. Given two different faces F\,F2, and an edge E, we calculate that 

' \£E^d(F 1 )nd(F 2 ), 
+1 if the "black to white" orientation of E 
Se-,Fi,F2 := { i s the same as the orientation induced (65) 

by the orientation of the face F\, 
— 1 otherwise. 

Taking now the sum over all edges E, we get 

£(Fi,F 2 ) :=e(d(F l ),d(F 2 )) =J2 S E;F lt F 2 . 

E 

Example. Given a face F and an oriented loop L on T, let us calculate e(F,L). Denote by 
Or£ the "black to white" orientation of the edge E. The orientation of L induces an orientation 
Or£ 1 of any edge E of L. Set 



Se,l = -Ot e ,l/Oi e e {±1}. 



(66) 
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So if the orientation of L is opposite to the canonical orientation of E, then 5e,l = 1- Otherwise 
5e,l = ~ 1- Then it follows from Figure [39] that we have 

e(F,L)= Yl 5 e,l. (67) 

EeLnd(F) 

Here the sum is over all edges E shared by the face F and the loop L. 

Changing the orientation of the surface amounts to changing the sign of the Poisson bracket. 
Interchanging black and white we change the sign of the Poisson bracket. 

8.3 Amalgamation and the Poisson structure on Cr- 

We extend the story to the case T is a bipartite graph on an oriented surface S with boundary. 
Pairs (r, S) are encoded by bipartite ribbon graphs. 

Bipartite ribbon graphs Take a bipartite graph embedded into an oriented surface S with 
boundary, so that some vertices of the graph may lie at the boundary of S. Let us add to every 
boundary vertex of the graph a few edges sticking out of the surface, called external legs. The 
graph r is then a bipartite ribbon graph. Every bipartite ribbon graph is obtained this way. 

Definition 8.4. Given a bipartite ribbon graph T, the moduli space C-p parametrizes line bundles 
with connection over the graph T, equipped with a trivialization of the restriction of the bundle 
to every external edge. 

Let Ti and T2 be bipartite ribbon graphs. Let us assume that we are given a way to glue 
some pairs of their external legs into a bipartite ribbon graph T. So each external leg of each 
graph is glued to no more than one external leg of the other graph, and the result is a bipartite 
ribbon graph. Then there is a map of moduli spaces, the amalgamation map: 

a : £ri x £t 2 — > £r- (68) 

We use trivializations at the external legs of the bundles on T\ and T2 to glue a new bundle on 
r. It inherits a connection. The following is obvious. 

Lemma 8.5. The gluing map I168\) is a principal fibration over Cy with the fiber (C*) , where 
k is the number of pairs of glued boundary legs. 





Figure 42: White and black hedgehogs with six external legs. 

A bipartite ribbon graph with single internal vertex is called a hedgehog, see Figure [421 It 
is a white / black hedgehog according to the color of its internal vertex. Let H v be a hedgehog 
with the internal vertex v. Then one has 

C Hv = Hom(A„,C*) = Ker(/u : (C*) m — > C*Y (69) 
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Here m is the number of legs and \x is the multiplication map. The first isomorphism assigns 
to an oriented path connecting two external legs the monodromy of the connection along this 
path. Since the line bundle is trivialized at the endpoints, it is just a non zero complex number. 
The second is obvious. 

A ribbon graph T is glued from the hedgehogs provided by the stars of its internal vertices. 
Amalgamating the hedgehog moduli spaces we get the moduli space Cr- 







Figure 43: Gluing bipartite ribbon graphs from hedgehogs. 

Thanks to identification (|69p . the skew symmetric local pairing on the abelian group A„ 
provided by the color of the internal vertex v of the hedgehog gives rise to a Poisson structure 
on the moduli space jCh v - The following proposition follows easily from the definitions. 

Lemma 8.6. There is a unique Poisson structure on the space Cy such that the Poisson structure 
on the space Ch v assigned to a hedgehog H v is the standard one, and the amalgamation maps 
Ii68\) are Poisson. 

When r is a bipartite graph on a surface without boundary we recover its Poisson structure. 
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